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A  COMPUTER  PROGRAM  FOR  ESTIMATION  OF  PARAMETERS  OF  THE  WEIBULL 
INTENSITY  FUNCTION  AND  FOR  THE  CRAMER-VON  MISES 
GOODNESS  OF  FIT  TEST 


1.  INTRODUCTION 


The  Weibull  intensity  function 


U(X)  =  X3X' 


(1.1) 


A>0,  B>0,  X>0,  is  frequently  used  as  a  model  for  the  determination  of 
reliability  growth  and  wear-out  characteristics  for  a  wide  variety  of 
complex,  repairable  systems.  The  failure  rates  of  military  equipment 
such  as  vehicles,  aircraft,  guided  missiles,  electronic  computer  systems, 
and  ammunition  are  being  evaluated  using  this  model. 

Formulas  have  been  developed  by  the  Army  Materiel  Systems  Analysis 
Activity  (AMSAA)  for  the  maximum  likelihood  estimation  of  the  unknown 
parameters  ^  and  6,  based  upon  sample  data.  These  estimation  formulas 
for  the  Weibull  process,  found  in  Crow  (1975),  can  be  stated  as 


X  ,nT?2lnT.2-T?,lnT„) 


K  N. 
i=l  j=l 


(1.2) 


where 


y  (T^  -T^  ) 

i2  'iV 


(1.3) 


K  is  the  number  of  systems  under  study; 

N.  is  the  total  number  of  failures  (or  occurrences  of  an  event 
'  under  study,  such  as  unscheduled  maintenance  actions;  etc.) 

'hh 

for  the  i  system; 

T-i  is  the  starting  time  of  the  period  of  continuous  observation 

I  i  u 

of  the  i^^  system; 

T.„  is  the  ending  time  of  the  period  of  continuous  observation 

1  ^  J.  L 

of  the  i^^  system; 

j.  u 

X. .  is  the  time  of  occurrence  of  the  failure  (or  event),  for 

'  J  th 

the  i  system; 


In  is  the  natural  logarithm,  and  O.ln(O)  is  defined  to  be  0. 
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To  expedite  the  computation  of  these  estimates,  AMSAA  developed  a 

^  As 

FORTRAN  computer  program  to  calculate  3  and  x.  That  program,  which  was 
documented  in  Belbot  (1974),  was  successfully  employed  by  the  U  S  Army 
Materiel  Development  and  Readiness  Command  (DARCOM),  various  subordinate 
commands  and  several  project  managers'  offices,  as  well  as  by  AMSAA. 

After  the  parameters  were  determined  by  that  program,  a  goodness  of  fit 
test  was  frequently  used  to  test  statistically  the  hypothesis  that  the 
failure  times  of  the  systems  being  analyzed  followed  a  nonhomogeneous 
Poisson  process  with  Weibull  intensity  function  (see  Crow  [1975]).  The 

modified  Cramer-Von  Mises  goodness  of  fit  statistic  was  computed,  either 
by  hand  or  by  a  separate  computer  program. 

Obviously,  the  consolidation  of  an  automated  goodness  of  fit  test  with 
the  computer  routine  which  estimates  the  parameters  3  and  x,  would  increase 
efficiency  and  accuracy.  Either  the  manual  calculations  or  the  use  of  a 
separate  program  would  be  eliminated.  To  accomplish  properly  this  consolida¬ 
tion,  the  estimation  procedure  was  subordinated  to  a  new  main  program  which 
also  controls  the  input  of  data  and  the  execution  of  the  goodness  of  fit 
test.  Because  of  the  radical  nature  of  this  redesign,  it  was  appropriate  to 
incorporate  other  new  features  at  the  same  time.  Principal  among  these  new 
features  are  a  simplified  input  procedure  and  dynamic  data  storage  alloca¬ 
tion.  The  resulting  computer  program  is  easier  to  use  and  provides  more 
information  that  its  antecedent  program.  This  note  will  explain  the 
structure  and  the  use  of  this  new  program. 

2.  COMPUTING  PROCEDURE 

2.1  Estimation  of  Parameters 


Since  the  formulas  (1.2)  and  (1.3)  do  not,  in  general,  yield 

A.  A 

3  and  x  in  closed  form,  an  iterative  technique  is  required.  Formula 
(1.2)  may  be  recast  as 


1\ 

I  I 

i=l  j=l 


K 


K  A 

I  (Tf2lnTi2-Tfiln(Tii) 
i=l _ 

1/  /V 

y(T^  -T^  ) 

^yi2  ‘ir 


(2.1) 


by  replacing  x  by  its  equivalent  expression  from  equation 
execution  of  a  few  simple  algebraic  operations.  Equation 

A 

consists  of  a  constant  with  regard  to  3,  minus  a  function 
0,  or  simply 

A  -  D(3)  =  0 


(1 .3) ,  and  by 
(2.1)  now 

A 

of  3,  yielding 
(2.2) 


The  correct  value  of  3  will  satisfy  equation  (2.2)  and  can  be  used  to 

A 

calculate  the  corresponding  value  of  x. 
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The  solution  for  3  is  iteratively  determined  in  the  following  way. 

For  an  initial  estimate  3  which  is  assumed  to  be  larger  than  the  true  3, 

A  I 

the  expression  A-D(3  )  is  evaluated.  For  all  values  of  3  larger  than 
the  true  3,  the  subtraction  yields  a  negative  result.  After  each  negative 

I 

result,  3  is  reduced  by  the  initial  step  size  of  1,  and  A-D(3  )  is  again 
evaluated. 

When  a  positive  number  results  from  the  subtraction,  indicating  that 

/s  f 

3  is  too  small,  the  step  size  is  decreased  to  0.10  of  the  present  step 

yv  I 

size,  the  previous  value  of  3  which  gave  a  negative  result  for 

A  I 

A-D(3  ),  is  adjusted  by  the  new  step  size  and  the  evaluation  process 
begins  again. 

A  I 

The  iteration  procedure  continues,  adjusting  3  by  the  new  step  sizes, 
until  the  left  side  of  equation  (2.2)  is  within  a  specified  tolerance  e  of 

0.  X  is  then  calculated  based  on  3,  using  factors  already  computed  in 

A 

finding  3.  This  procedure  is  summarized  by  the  state  diagram  (Figure  2.1). 
2.2  Goodness  of  Fit  Test 

The  Cramer-Von  Mises  Goodness  of  Fit  Test  is  appropriate  whenever 
the  starting  time  of  each  system  is  equal  to  0.  To  perform  this  test, 
the  program  first  transforms  the  failure  times.  For  time  truncated  test¬ 
ing,  the  failure  times  for  each  system  are  divided  by  the  ending  time  of 
the  test  period  for  that  system.  In  failure  truncated  testing,  for  every 
system,  all  the  failure  times  except  the  last,  are  divided  by  the  last 
failure  time.  The  last  failure  time  is  thereafter  excluded  from  the 
calculations  and  the  number  of  transformed  failures  is  one  less  than  the 
original  number  of  failures  for  each  such  system.  All  the  transformed 
failure  times  are  then  sorted  into  increasing  order. 

^  Next,  the  unbiased  estimate  3  of  the  estimated  shape  parameter 

3,  is  calculated  using  the  relation: 

B  =  ^  3  (2.3) 


where: 

M  is  the  number  of  transformed  failure  times; 
and  N  is  the  number  of  original  failure  times. 
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ENTER 


EVALUATE  A 


Figure  2. 1  State  Diagram 
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formula: 


Finally,  the  Cramer-Von  Mises  statistic  C 


2 

M 


is  computed  by  the 


1 


M  12M 


M 

I 

i=l 


(Z^  -  ilzi) 
^^i  2M  ’ 


(2.4) 


where  the  Z.  are  the  transformed  failure  times.  An  explanation  of  this 

1  ip 

test  and  a  table  of  critical  values  of  Cj^  may  be  found  in  Crow  (1975). 


Next  page  is  blank. 


5 


3.  DESCRIPTION  OF  PROGRAM 


3.1  Major  Features 

The  program,  which  is  listed  in  Appendix  A,  has  some  important 
features.  First,  the  program  is  written  in  American  National  Standards 
Institute  (ANSI)  FORTRAN  X3. 9-1966,  and  should  therefore  execute  on  any 
computer  having  a  compiler  for  this  language.  Secondly,  the  amount  of 
storage  required  to  use  this  program  should  not  cause  difficulties  since 
all  data  arrays  are  dynamically  allocated  under  control  of  the  main 
program  (see  Chung-Phil lips,  et  al.,  [1975]). 

Finally,  and  perhaps  most  importantly  for  the  user,  this  program 
uses  free-field  input,  that  is,  no  specific  format  is  required  for  the 
input  data.  The  time  and  the  effort  regularly  expended  in  preparing 
data  for  input,  are  greatly  reduced  because  of  this  feature.  Moreover, 
the  misalignment  of  data  fields  to  formats,  a  frequent  source  of  errors 
in  using  many  computer  programs,  is  eliminated  entirely.  While  no  input 
scheme  can  be  regarded  as  foolproof,  free-field  input  is  much  more  flexible 
than  fixed-field  format  specifications. 

3.2  Overall  Characteristics 


All  calculations  in  the  program  are  made  in  double  precision  mode. 
Experience  has  shown  that  the  use  of  single  precision  variables  for  these 
calculations  often  results  in  significant  discrepancies  in  the  estimates  of 
the  parameters  due  to  errors  accumulated  during  the  iterative  process. 

The  modular  structure  of  the  program  (see  Figure  3.1)  reflects 
organization  by  functional  purpose.  The  input  of  data,  certain  inter¬ 
mediate  calculations,  the  estimation  of  parameters,  and  the  goodness  of 
fit  test  are  each  performed  by  an  independent  module.  Major  subprograms 
print  their  results  as  the  values  become  available.  Subroutines  which 
detect  errors,  print  diagnostic  messages  naming  the  detecting  routine 
and  briefly  stating  the  difficulty,  and  then  attempt  to  continue  process¬ 
ing  when  possible.  Independence  of  the  subroutines  is  maintained  by 
restricting  communication  between  individual  subprograms  to  the  passage 
of  formal  parameters  in  argument  lists.  No  COMMON  statements  are  used. 

3.3  Specific  Details  of  Routines 

In  addition  to  controlling  the  major  modules,  the  main  program 
also  allocates  storage  for  data  arrays,  as  stated  in  Section  3.1.  The 
allocation  is  based  upon  the  maximum  number  of  failures,  NFAIL,  and  the 
maximum  number  of  systems,  NSYS.  The  master  data  array,  BLOCK,  has  a 
length  of  NTOTAL,  equal  to  the  value  of  NFAIL  plus  six  times  the  value 
of  NSYS.  To  redimension  the  entire  program,  one  merely  adjusts  the 
values  of  NFAIL,  NSYS  and  NTOTAL  in  the  DATA  statements  at  the  beginning 
of  the  main  routine  and  changes  the  size  of  the  array  BLOCK,  also  found 
at  the  beginning  of  the  main  program,  to  equal  the  new  value  of  NTOTAL. 
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INDATA  LOGS  MLEW  CVM 


Figure  3.1  Module  Organization  Chart. 


By  these  actions,  all  data  arrays  in  all  subroutines  will  be  properly 
resized.  Since  the  program  size  excluding  data  arrays,  is  less  than^ 

5,000  words,  computer  memory  requirements  can  be  scaled  to  problem  size 
through  use  of  this  feature. 

Also  found  in  the  DATA  statements  at  the  beginning  of  the  main 
routine  are  the  unit  number  for  input,  lUNIT,  and  the  logical  switch  ECHO 
which  controls  the  printing  of  the  input  data.  Just  as  the  storage 
allocation  values,  these  values  may  be  changed  as  needed. 

The  first  major  module  is  for  the  input  of  data.  The  INDATA 
module  reads  the  beginning  and  the  ending  times  and  the  failure  times 
from  the  input  unit,  lUNIT.  If  the  logical  variable  ECHO  is  true,  the 
submodule  REPEAT  will  print  the  input,  using  the  LOCIP  subroutine  to 
isolate  in  storage  the  failures  for  each  system  and  the  PRNTX  subroutine 
to  print  them.  The  failure  times  are  stored  in  a  linear  array  with  nega¬ 
tive  signs  appended  to  the  failures  associated  with  even  numbered  systems. 
This  scheme  preserves  the  identification  of  the  failures  with  the 
respective  systems  without  using  additional  storage. 

The  free-field  reading  of  data  is  performed  by  the  INREAL  sub- 
module.  This  submodule,  consisting  of  the  routines  INREAL,  CHAR,  INSYML, 
MATCH  and  UNSGND,  is  a  translation  from  ALGOL  into  FORTRAN  of  Algorithm 
239  of  the  Association  for  Computing  Machinery  (see  McKeeman  [1964]). 
Specific  details  concerning  the  input  arrangement  are  given  in  Section  4. 

The  INDATA  subroutine  also  sets  three  logical  variables  depend¬ 
ing  on  the  input.  If  data  errors  are  encountered,  the  variable  NOGOOD  is 
made  true.  If  all  systems  start  at  0.0,  then  the  goodness  of  fit  test 
will  be  appropriate  and  so  the  variable  GOF  is  set  to  true.  Lastly,  when 
the  end  of  input  is  reached,  the  logical  variable  HALT  is  returned  as  true 

The  second  principal  module,  LOGS,  calculates  the  logarithms  of 
the  beginning  and  the  ending  times,  and  the  sum  of  the  logarithms  of  the 
failure  times.  For  computational  purposes,  a  beginning  time  of  0.0  is 
defined  to  have  a  logarithm  of  0.0  instead  of  infinity.  Note  that 
failure  times  of  0.0  are  not  valid  for  this  model  and  are  flagged  as 
errors  by  the  preceding  INDATA  module. 

The  next  major  module,  MLEW,  computes  the  maximum  likelihood 
estimates  of  the  parameters  of  the  Weibull  intensity  function  for  the 
given  data,  using  the  formulas  discussed  in  Section  2.1.  If  unsuccessful, 
this  module  will  report  one  of  three  possible  error  conditions.  The 
first  error  message  "BETA  LESS  THAN  0.00000001"  indicates  that  the 

A 

data  should  be  rechecked  for  the  reasonableness  of  a  very  small  f.  The 
second  error  message  is  "INITIAL  ESTIMATE  OF  BETA  IS  TOO  SMALL."  Since 
the  initial  estimate  of  BETA  is  set  to  10  at  the  start  of  the  MLEW  sub¬ 
program,  this  message  indicates  some  pecularity  of  the  data.  (In 

general,  0<3<10.)  The  third  message,  "STEP-SIZE  HAS  BECOME  INSIGNIFICANT 
BETA  NOT  RESOLVABLE,"  indicates  that  the  module  has  gone  as  far  as 
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possible  trying  to  meet  the  tolerance  set  for  the  difference  A-D(6  ). 

This  tolerance,  EPSILN,  may  be  enlarged  by  changing  the  assignment  state¬ 
ment  also  located  at  the  beginning  of  the  MLEW  subprogram. 

The  last  major  module,  CVM,  performs  the  Cramer-Von  Mises  goodness 
of  fit  test,  as  described  in  Section  2.2.  The  failure  times  for  each 
system  are  located  in  storage  using  the  LOCIP  subroutine,  and  examined  to 
determine  if  the  testing  was  time  truncated  or  failure  truncated.  The 
failure  times  are  then  transformed  by  the  TRANSF  subroutine  and  sorted  by 
the  SORTX  subroutine.  (SORTX  is  a  modification  of  an  utility  subprogram 
described  in  Campbell,  et  al.,  [1970].)  The  unbiased  estimate  of  3,  UNBETA, 
is  calculated  next,  as  explained  earlier.  The  last  phase  depends  on  the 
system  starting  times.  If  any  starting  time  is  non-zero,  the  module 

terminates  with  a  message  stating  that  the  Cramer-Von  Mises  goodness  of 
fit  test  is  not  appropriate.  Otherwise,  the  goodness  of  fit  statistic, 

CM2,  is  computed  and  printed. 
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4.  INPUT  REQUIREMENTS 


As  stated  previously,  the  input  for  this  program  is  free-field.  The 
only  requirement  regarding  spacing  is  that  at  least  one  blank  column 
separate  adjacent  values.  The  values  must  not  run  together.  This  means 
that  the  program  generally  takes  the  same  view  of  the  data  that  a  person 
would,  namely,  that  each  cluster  of  numeric  characters  constitutes  one 
data  value.  The  only  exception  to  this  rule  occurs  at  the  boundaries  of 
records.  Since  the  input  is  treated  as  a  continuous  stream,  a  string  of 
characters  beginning  in  the  first  column  of  a  record,  is  considered  a 
continuation  of  the  string  of  characters  ending  in  the  last  column  of  the 
previous  record,  if  any.  Record  boundaries  are  not  delimiters;  blanks 
are  the  only  delimiters. 

The  data  required  for  this  program  consist  of  the  beginning  and  ending 
times  of  the  test  period  for  each  system,  and  the  failure  times  for  each 
system.  The  arrangement  of  the  input,  (which  is  also  stated  in  the 
comments  of  the  INDATA  subroutine),  is  by  system.  The  first  data  value 
is  the  beginning  time  of  the  first  system.  The  second  data  value  is 
the  ending  time  of  the  first  system.  Next  is  the  failure  times  for  the 
first  system,  followed  by  a  negative  value  to  mark  the  end  of  the  first 
system.  The  same  pattern,  beginning  time,  ending  time,  failure  times, 
and  negative  trailer,  is  repeated  for  each  subsequent  system  in  the  first 
data  case. 

Another  negative  value  (making  two  in  a  row),  signals  the  end  of 
input  of  the  current  data  case,  and  the  beginning  of  the  computational 
procedures.  The  same  arrangement  may  be  repeated  for  as  many  cases  as 
desired  per  program  run.  When  the  input  routine  encounters  a  negative 
value  after  completing  a  case,  (that  is,  the  third  negative  value  in  a 
row),  the  end  of  the  program  run  is  indicated. 

Thus,  as  a  simple  example,  if  one  desired  to  use  this  program  for 
one  run  consisting  of  one  case  wherein  one  system  experienced  seven 
failures,  the  input  data  would  be:  the  beginning  time,  the  ending 
time,  the  seven  failure  times,  and  three  negative  values.  To  demonstrate 
the  latitude  of  the  input  requirements,  the  data  for  a  number  of  test 
cases  are  shown  in  Appendix  B.  Notice  that  any  negative  value  is  accept¬ 
able  as  a  trailer  and  that  data  values  may  be  entered  with  or  without 
decimal  points.  Although  not  shown  in  the  examples,  data  values  may 
also  be  in  exponential  form,  that  is,  containing  'E,'  '+'  or 
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5.  TRANSFERABILITY  AND  MODIFICATION 


Since  this  program  is  written  in  standard  FORTRAN,  transfer  to  other 
computer  systems  should  be  straightforward.  To  assist  in  the  transfer 
process.  Appendix  C  contains  the  output  produced  by  the  program  for  the 
input  shown  in  Appendix  B.  This  output  was  generated  on  a  Control  Data 
Corporation  (CDC)  Cyber  76  Computer,  using  the  program  exactly  as  listed 
in  Appendix  A.  (Note  that  non-standard  PROGRAM  statement  required  by 
the  CDC  Cyber.) 

The  input  for  these  test  cases  came  from  records  of  eighty  characters 
each.  If  the  input  record  length  is  other  than  eighty,  two  changes  may 
be  required.  The  value  of  the  variable  LENGTH  and,  if  necessary,  the 
dimensioned  size  of  the  array  BUFFER,  should  be  adjusted  in  the  INSYML 
subroutine. 

Alternatively,  one  could  replace  the  entire  INREAL  submodule.  Although 
these  routines  were  written  to  be  fully  transportable,  running  time  might 
be  saved  by  using  the  system  defined  free-field  reading  capability  of  any 
computer  having  such  a  feature.  As  an  example,  a  substitute  for  the 
INREAL  submodule,  suitable  for  the  CDC  Cyber  76,  is  shown  in  Figure  5.1. 
Such  substitutes,  however,  are  system  dependent  and  not  readily  trans¬ 
ferable. 
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FIGURE  5.1  SUBSTITUTE  FOR  INREAL  SUBMODULE 


SUBROUTINE  INREAL  (lU,  X) 

DOUBLE  PRECISION  X 

READ  (lU,  *)  X 

RETURN 

END 
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LISTING  OF  PROGRAM 
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PPOGPAH  MAIN  (INPtT-/8CJ»0UTPUT»TAPE5-INPUT»TAPE6-0UTPUT)  MAIN  10 

THE  ABOVE  STATEMENT  IS  NON-STANDARD#  BUT  REQUIRED  FOR  MAIN  20 

COC  FORTRAN,  MAIN  30 

MAIN  40 
MAIN  50 

MAIN  DRIVER  FOR  UEIBULL  INTENSITY  MODEL  PARAMETER  MAIN  60 

ESTIMATION  AND  GOODNESS  OF  FIT  TEST.  MAIN  70 

MAIN  eo 

(VERSION  OF  27  JULY  197Q)  MAIN  90 

MAIN  100 

INPUT  REQUIREMENTS  ARE  DESCRIBED  IN  THE  'INDATA*  MAIN  110 

SUBROUTINE.  MAIN  12C 

MAIN  130 
MAIN  140 

DOUBLE  PRECISION  BETA#  SUMLNX  MAIN  150 

DOUBLE  PRECISION  BLOCKdlOOO)  MAIN  160 

LOGICAL  ECHO#  FAULT#  GCF#  HALT#  NOGOOO  MAIN  170 

MAIN  180 

DATA  lUNIT  /5/  main  190 

DATA  ECHO  /.TRUE./  MAIN  20C 

DATA  NEAR#  NSYS  /530C#  1000/  MAIN  21C 

DATA  NTOTAL  /lltyv./  MAIN  220 

MAIN  230 

WRITE  (6»20)  main  240 

MAIN  250 

ALLOCATE  STORAGE  BASED  ON  MAXIMUM  NUMBERS  OF  FAILURES  MAIN  260 

AND  SYSTEMS.  MAIN  270 

NS2«NSYS*2  main  280 

Il«l  MAIN  290 

12-11+NFAIL  MAIN  300 

I3«I2+NS2  main  310 

I4-I3+NS2  main  32C 

IT0TAL-I4+NS2-1  MAIN  330 

IF  (ITOTAL.LE. NTOTAL)  GO  TO  K  MAIN  340 

WRITE  (6#30)  ITGTAL/NTCTAL  MAIN  350 

STOP  main  360 

MAIN  370 

BEGIN  PROCESSING.  MAIN  380 

1C  CALL  INOATA  ( BLOCK ( 11 )» BLOCK ( 12 )# NOGOOD# GOF# HALT.NF AIL# NSYS » M#K# lUM A  IN  390 

lNn#ECHn)  MAIN  40C 

IF  (HALT)  STOP  MAIN  410 

If  (NDGOOD)  GO  TO  1C.  MAIN  420 

CALL  LOGS  (PL0CK(I1)#SL0CK(I2)#9L0CK(I4)#NFAIL#NSYS#M#K#SUMLNX)  MAIN  430 

CALL  MLEW  (SUMLNX#BLCCK(I2)#PL0CK(I3)#BL0CK(I4)»FAULT#NSYS#M,K#PETMAIN  44C 

lA)  MAIN  450 

IF  (.NOT. FAULT)  CALL  CVM  ( BL OCK ( 1 1 ) » B LOCK ( 1 2 ) , NF A I L » NS YS» M, K# 6E T A# M A I N  460 

IGOF)  MAIN  47C 

GO  TO  10  MAIN  480 

main  490 

20  FORMAT  (71H1  W'EIBULL  INTENSITY  MODEL  PARAMETER  ESTIMATION  AND  GOODMAIN  500 

INESS  OF  FIT  TEST/25HC  VtRSION  OF  27  JULY  1979////)  MAIN  510 

30  FORMAT  (33H1  AMOUNT  OF  STORAGE  RECUESTED  IS  #I6»40H  WORDS.  AMOUNTMAIN  520 

1  OF  STORAGE  AVAILABLE  IS  #I5#25H  WORDS.  PROGRAM  ABORTED,/)  MAIN  530 

END  main  540 


MAIN  DRIVER  FOR  WEIBULL  INTENSITY  MODEL  PARAMETER 
ESTIMATION  AND  GOODNESS  OF  FIT  TEST. 

(VERSION  GE  27  JULY  1979) 

INPUT  REQUIREMENTS  ARE  DESCRIBED  IN  THE  'INDATA* 
SUBROUTINE. 


DOUBLE  PRECISION  BETA#  SUMLNX 
DOUBLE  PRECISION  BLOCKdlOOO) 

LOGICAL  ECHO#  FAULT#  GCF#  HALT#  NOGOOO 

DATA  lUNIT  /5/ 

DATA  ECHO  /.TRUE./ 

DATA  NFAIL#  NSYS  /5D0C#  1000/ 

DATA  NTOTAL  /lltyv./ 

WRITE  (6»20) 

ALLOCATE  STORAGE  BASED  ON  MAXIMUM  NUMBERS  OF  FAILURES 
AND  SYSTEMS, 

NS2«NSYS*2 

11- l 

12- I1+NFA1L 

13- I2+NS2 

14- I3+NS2 
IT0TAL-I4+NS2-1 

IF  (ITOTAL.LE. NTOTAL)  GO  TO  K 
WRITE  (6#30)  ITGTA L#NTCTAL 
STOP 
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SUBROUTINE  CVM  ( X*C»NFAIL»NSYS»M»K>BETA»GOF ) 

COMMENT  THIS  SUBROUTINE  PERFORMS  THE  CRAMER-VDN  MISES  GOODNESS 

C  UF  FIT  TEST. 

C  f 

OCUBLE  PRECISION  X(NFAIL»»  C(NSYS,2I 

DOUBLE  PRECISION  BETA,  CM2»  OM»  SUMSQS»  TO»  TWOM»  UNBETA 
DOUBLE  PRECISION  TERM,  TERMl,  TEBM2 
LOGICAL  GOF 
LOGICAL  TINETR 
C 

WRITE  (6,bO)  BETA,M 

IP«G 

IBT-1 

TIHETR«,TRUE. 

N«M 

C 

DO  30  J»1»K 
I3«IP*1 

CALL  LOCIP  (X,NFAIL,N,IP> 

IE«IP 

DO  1C  I«IB»IE 

C  CHECK  FOR  FAILURE  TRUNCATED  TESTING. 

IF  (DABS (DABS! X{ 1) )-C ( J, 2 ) I .LE.l.t 0-081  TIMETR* . FALSE . 

10  CONTINUE 

IF  (TIMETR)  GO  TO  20 

IE-IE-1 

M»  M—  1 

TIHeTR«.TRUE . 

20  CONTINUE 
TO»C( J, 2) 
iET«ie 

C  TRANSFORM  THE  FAIIURES. 

CALL  TkANSF  ( X» TO, IB, IBT,  IET.NFAIL ) 

IBT-IET+l 
30  CONTINUE 
C 

DO  4'J  I»1,M 
X(I)*OABS(X(  D) 

4C  CONTINUE 

C  SORT  THE  TRANSFORMED  FAILURES  INTO  INCREASING  ORDER. 

CALL  SORTX  (X,M) 

C 

DM*D6LE(FLQAT(M)  ) 

TwOM*2.CDC*OM 

C  UNBIASED  estimate  OF  BETA. 

UN3ETA-BETA* (DM-l.OOC )/OBLE(FLOAT(N)) 

WRITE  (6,11D)  UNBETA 
IE  (COF)  GO  TO  50 
WRITE  (6,133) 

GO  TO  70 
C 

6't  SUMSOS-C.CDO 
DO  60  I»1,M 
TEPM1»X(I )**UNBETA 
TtRM2«DBL£(FLQAT(2*I-l))/TW0M 
TERM.TERM1-TERM2 
SUMS0S«SUMSQS*T£RM**2 


CVM 

10 

CVM 

20 

CVM 

30 

CVM 

AO 

CVM 

5C 

CVM 

6C 

CVM 

7C 

CVM 

80 

CVM 

90 

CVM 

100 

CVM 

110 

CVM 

120 

CVM 

130 

CVM 

lAO 

CVM 

15C 

CVM 

160 

CVM 

17C 

CVM 

180 

CVM 

19C 

CVM 

200 

CVM 

210 

CVM 

22C 

CVM 

230 

CVM 

2A0 

CVM 

250 

CVM 

26C 

CVM 

270 

CVM 

280 

CVM 

290 

CVM 

300 

CVM 

3  1C 

CVM 

320 

CVM 

330 

CVM 

340 

CVM 

350 

CVM 

360 

CVM 

3  70 

CVM 

380 

CVM 

390 

CVM 

400 

CVM 

41C 

CVM 

420 

CVM 

430 

CVM 

440 

CVM 

450 

CVM 

460 

CVM 

470 

CVM 

460 

CVM 

490 

CVM 

500 

CVM 

510 

CVM 

520 

CVM 

530 

CVM 

540 

CVM 

550 

CVM 

560 

CVM 

5  7C 

22 


60 


CONTINUE 

CRAMER-VON  MISES  STATISTIC. 
CN2-SUHSQS+(1.0D0M12.C0C*DM)  ) 
WRITE  C6#90>  CM2#M 


70  RETURN 


ao 


CVf^ 
CVM 
CVM 
CVM 
CVM 
CVM 
CVM 

FORMAT  <42H1  CRAMER  -  VON  MISES  GOODNESS  OF  FIT  TEST. // /18H0EST IMACVM 
ITEP  BETA  •  #iPDli).7/22H0NUM6ER  OF  FAILLRES  ■  #16)  CVM 

90  FORMAT  (32H0CRAMER  -  VON  MISES  STATISTIC  •  # 1P015 . 7/52H0R EJECT  THECVM 
1  WEIBULL  INTENSITY  MODEL  IF  THE  STAT ISTIC /48H  EXCEEDS  THE  APFROPRICVM 
2ATE  CRITICAL  VALUE  FOR  M  •  #U/)  CVM 

1-JO  FORMAT  (77HOTHE  CRAMER  -  VON  MISES  GOODNESS  OF  FIT  TEST  IS  NOT  APPCVM 
IROPPIATE  FOR  THIS  CASE/58H  BECAUSE  ONE  OR  MORE  SYSTEMS  HAVE  NON-ZECVM 
2R0  STARTING  TIMES./)  CVM 

110  format  (29H0UNBIASE0  ESTIMATE  OF  BETA  •  f 1PD15.7//// )  CVM 

END  CVM 


580 

590 

600 

610 

620 

630 

640 

650 

660 

670 

680 

690 

700 

710 

720 

730 

74C 
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SUBROUTINE  INOATA  (X>C*N06000»GOF»HALT>NFAIL»NSYS»M»K»IUN1T»ECHO) 
COMMENT  THIS  SUBROUTINE  READS  IN  THE  BEGINNING  AND  ENDING  TIMES» 

C  AND  THE  FAILURE  TIMES. 

THE  ARRANGEMENT  OF  INPUT  IS  AS  FOLLOWSI 

BEGINNING  AND  ENDING  TIMES  FOR  FIRST  SYSTEM.  FAILURE 

TIMES  FOR  FIRST  SYSTEM.  NEGATIVE  VALUE  AS  TRAILER. 
BEGINNING  AND  ENDING  TIMES  FOR  SECOND  SYSTEM.  FAILURE 

TIMES  FOR  SECOND  SYSTEM.  NEGATIVE  VALUE  FOR  TRAILER. 


«  «  • 

•  •  • 

BEGINNING  AND  ENDING  TIMES  FOR  K-TH  SYSTEM.  FAILURE 

TIMES  FOR  K-TH  SYSTEM.  NEGATIVE  VALUE  AS  TRAILER. 
negative  VALUE  TO  MARK  END  OF  CASE. 

(REPEAT  ABOVE  FOR  AS  MANY  CASES  AS  NEEDED.) 

negative  value  to  mark  end  of  run. 

INPUT  IS  FREF-FIELD.  REQUIRING  ONLY  THAT  AT  LEAST  ONE  BLANK 
COLUMN  SEPARATE  ADJACENT  VALUES. 


C 


C 


C 

C 

c 

c 


c 

c 


c 

c 


DOUBLE  PRECISION  X(NFAIL).  C(N5YS.2) 
LOGICAL  ECHO.  GOF.  HALT.  NOGOOD 

GCF».TPUE. 

HALT-. FALSE. 

NOGOOD- . FALSE . 

WRITE  (6.130) 


J-1 

I-l 


BEGIN  INPUT  CYCLE. 


10  CALL  INPEAL  (  I UM  T  »  C  (  J  .  1 )  ) 

NEGATIVE  VALUE  TO  MARK  THE  END  OF  THIS  CASE. 

IF  (C(  J.D.LT.C.CDO)  60  TO  50 
CALL  INPEAL  (IUNIT.C(J.2)) 

IF  (OABS(C( J.i) ) .GT.1.0D-C8)  GOF-. FALSE. 

J*  J  +  1 

IF  ( J.Lt.NSYS)  GO  TO  20 
WRITE  (6.100)  J 
NOGOOD- . TRUE . 

20  CALL  IMREAL  (  IUNIT.X( I ) ) 

NEGATIVE  VALUE  TO  MARK  THE  END  OF  THIS  SYSTEM. 

IF  (X( I  )  .LT.0.000)  CO  TO  10 
IF  (X(  I  )  .GT. 1.00-15)  GO  TO  30 
WRITE  (6.110) 

NOGOOD-. TRUE. 

EACH  FAILURE  MUST  FALL  WITHIN  THE  TEST  pEPIOO. 
30  IF  (C ( J-1. 1 ) . LE . X( I ) . ANO.X ( I ) . L6 .C ( J-1. 2 ) )  GO  TO  40 
WRITE  (6.120)  X ( I ) .C ( J-1. 1 )»C ( J-1. 2 ) 

NOGOOD- .TRUE . 

4C  IF  (MOO(J.2).EC.O)  X(I)»-X(I) 


INOAA  10 

INDAA  2C 

INOAA  30 

INDAA  40 

INDAA  50 

INOAA  60 

INOAA  70 

INOAA  SO 

INOAA  9C 

INOAAIOC 

INDAAllC 

INDAA12C 

INDAA130 

INDAA140 

IN0AA150 

IN0AAi60 

INOAA170 

INDAA160 

IN0AA190 

INDAA200 

INDAA21U 

INOAA22C 

INDAA230 

IN0AA24O 

INDAA250 

IN0AA260 

INDAA27C 

INDAAZeC 

IN0AA290 

INOAA300 

IN0AA31C 

INDAA320 

IN0AA33O 

INDAA34C 

IN0AA35C' 

IN0AA360 

IN0AA370 

IN0AA38;', 

IN0AA39;< 

INOAA400 

IN0AA41C 

INDAA420 

IN0AA430 

INDAA440 

IN0AA45') 

INDAA460 

IN0AA47 J 

INDAA4EC 

IN0AA490 

INDAA50C 

IN0AA510 

IN0AA520 

INDAA530 

INDAA540 

IN0AA55C 

INDAA560 

IN0AA570 
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IN0AA580 

IF  (I.LE.MFAILI  GO  TC  20  IN0AA59C 

WRITE  I6»90»  I  INDAA6C0 

NOGOllO-.TRUE.  INDAA610 

GO  TO  2C  IN0AA62C 

C  -  IN0AA630 

C  END  OF  INPUT  CYCLE.  IN0AA640 

- - INDAA65J 

M* I— I  INDAA660 

IF  (P.LE.O)  GO  TO  6C  IN0AA670 

C  TOTAL  NUMBER  OF  FAILURES.  IN0AA68C 

WRITE  16,70»  M  INDAA690 

K«J>I  INDAA70C 

C  TOTAL  NUMBER  OF  SYSTEMS.  IN0AA710 

WRITE  (6»60I  K  IF:DAA720 

IF  (bCriOI  CALL  REPEAT  IX,C»NFAIL»NSYS»M,K)  IN0AA73C 

RETLRN  ■  IN0AA7A0 

- IN0AA75?' 

END  OF  FUN.  IN0AA760 

-  INDAA77C 

fcO  HALT«.IFUE.  INDAATBO 

WRITE  (6»140»  INDAA790 

return  IN0AA8C0 

C  INOAASiC 

r:  FORMAT  (ZeHCTOTAL  NUFflER  OF  FAILURES  •  »I5»  INDAA820 

bJ  FORMAT  (28HCTOTAL  NUMBER  OF  SYSTEMS  ■  »I5)  INDAABBC 

90  FORMAT  (18H0  INOATA  ERROR  —  ,I6,22H  IS  TOO  MANY  FAILURES./)  INOAA04) 

100  FORMAT  (IbHO  INDATA  ERROR  —  »I6»21H  IS  TOO  MANY  SYSTEMS./)  lNDAA85vT 

IIG  FORMAT  <52H0  INOATA  ERROR  —  A  FAILURE  AT  0.0000000  WAS  INPUT.  /71INDAA86U 
14C  The  PROBABILITY  OF  SUCH  A  FAILURE  TIME  IS  0.0  ACCORDING  TO  THE  IN0AAB7C 

ZMODEL./)  INOAASeO 

120  FORMAT  <33HC  INOATA  ERROR  —  THE  FAILURE  AT  »1PD10.3/44H  DOES  NOTINDAAeQO 

1  FALL  WITHIN  THE  TEST  PERIOD  FROM  ,1PD10.3»4H  TO  .1P010.3/)  IN0AA900 

130  FERMAT  (19H1  DATA  INPUT  PHASE.//)  INDAA91C 

140  FORMAT  (28HU  PROGRAM  RUN  ENOS  NORM ALL Y . / 1  HI )  INDAA920 

end  INDAA930 
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SUftRQUTINt  LOCIP  ( )(»NFAIL»M»  IP  » 


LQCIP  10 


CUhMLNT 

C 


THIS  SUBROUTINb  LOCATES  THE  POSITION  OF  THE  LAST  FAILURE  LOCIP  20 


associated  with  the  system  which  had  FAILURE  'XTIPJ*. 
RETURN»  'IP*  INDEXES  THIS  LAST  FAILURE. 


DOUBLE  PRECISION  X<NFAIL> 


LOCIP  30 
LOCIP  AO 
LOCIP  50 
LOCIP  60 
LOCIP  7.1 
LOCIP  8C 
LOCIP  90 
LOCIPlOO 


IB-lP  +  1 

LOCIP  90 

DoSo  1-IB,MH1  LOCIPlOO 

IF  ((x(I).LT.O.OOO,ANO.XH+1>.LT.O.OOO).OR.(X(I).GE.O.COO.ANO.X(I*LOCIP11C 

ID.GE.O.DOO)  )  GO  TO  10  L0CIP130 

pHuRN  LOCIPIAO 

1C  CONTINUE  tSciPlOO 

^  LOC IP170 

[SciPlBO 

t  Nl? 


% 
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SUBROUTINE  LOGS  ( X,Cf CLN, NF AIL#NSYS#M,K,SUMLNX ) 

CUMMENT  THIS  SUBROUTINE  TAKES  LOGARITHMS. 

C 

DOUBLE  PRECISION  X(NFAIL)>  C(NSYS#2)f  CLN(NSYS#2I 
DOUBLE  PRECISION  SUMLNX 
C 

C  LOGARITHMS  OF  BEGINNING  AND  ENDING  TIMES, 

DO  20  I«1,K 
CLN(I,2)«DL0G<C<I.2)) 

IE  (C(  If  D.LE.O.ODO)  GO  TO  10 
CLNdf  l)*DLQG(C<If  in 
GO  TO  2f 
1')  CUfD-O.OOO 
CLNdf  i)«0,0D0 
2'^  CONTINUE 
C 

C  SUM  OF  LOGARITHMS  OF  FAILURES, 

SUMLNX*O.COO 
DC  30  l«lfM 

SUMLNX«SUMLNXtDLQG(DA8S  (X(  Id) 

31  CONTINUE 
RETURN 
C 

END 


LOGS  1C 
LOGS  2C 
LOGS  30 
LOGS  40 
LOGS  5v/ 
LOGS  6C 
LOGS  7D 
LOGS  BO 
LOGS  ^0 
LOGS  iOu 
LOGS  110 
LOGS  120 
LOGS  130 
LOGS  140 
LOGS  ISO 
LOGS  It: 
LOGS  17C 
LOGS  IBO 
LOGS 
LOGS 

LOGS  210 
LOGS  22C 
LOGS  230 
LOGS  2^^ 
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SUBROUTINE  HLEW  ( SUMLNX# C» CB» CLN» N0G000»NS YS»H»K»BETA  ) 


ESTIMATES  OF  PARAMETERS  OF  THE  WEIBULL  INTENSITY  MODEL. 


THIS  SUBROUTINE  ESTIMATES  BETA  AND  LAMBDA  OF  THE 
WEIBULL  INTENSITY  FUNCTION  R (X» -LAMBDA*BETA*X** (BETA-1. C » 

BETA  IS  DETERMINED  BY  AN  ITERATIVE  PROCESS  WHICH 
EXAMINES  THE  SIGNED  DIFFERENCE  OF  A  CONSTANT  MINUS  A  FUNCTION  OF 
BETA  AS  BETA  IS  DECREASED  FROM  A  LARGE  INITIAL  ESTIMATE  BY 
NON-POSITIWE  POWERS  OF  10.0  UNTIL  AN  EPSILON  TOLERANCE  IS 
SATISFIED  OR  BETA  IS  LESS  THAN  0.00000001  IN  VALUE. 

LAMBDA  IS  CALCULATED  BASED  ON  BETA. 


DOUBLE  PRECISION 
DOUBLE  PRECISION 
DOUBLE  PRECISION 
DOUBLE  PRECISION 
DOUBLE  PRECISION 
LOGICAL  NOGOOD 


C(NSYS»2)»  CB(NSYS»2)# 
BETA#  LAMBDA#  EPSILN 
OENOM#  DENOMl#  OENOM2# 
A#  ABDI#  AOJ#  D#  DIFF# 
SUMLNX 


CLN(NSYS#2) 

TOP#  TOPI#  T0P2 
TOTFAL 


WRITE  (6#11C) 

BETA-i. 00*01 

EPSILN-l.OD-05 

TPTFAL«DBLE<FLOAT{M)) 

NPDIFF-0 

NOGOOD-. FALSE. 


CONSTANT  NOT  INVOLVING  BETA. 
A-SL’MLNY/TOTFAL 
WRITE  (5#130)  A 


BEGIN  ITERATION  PROCEDURE. 


ADJ-1.000 
IZEPO'O 
WRITE  (6#15C) 

C 

10  TDFl-C.ODO 
TCP2*0.00C 
DEMOMl-C.ODC 
penomb-o.odo 

DO  20  I«1#K 
DO  20  J-l#2 
CB(I#J)-C(I#J)**BETA 

20  CONTINUE 
C 

DO  BO  I-1#K 

TOPl-TOPl*CB(I#l)*CLN{I»i) 
T0P2«TnP2*CB(I#2)*CLN(I#2» 
DEKOMl-OENOMl  +  CBd#!) 
DENOM2«0£NOM2*CB(I#2) 

30  CONTINUE 
C 

TDP-T0P2-TCP1 

DEN0M-DENOM2-0£NaMl 

D« (TCP/DENCM>-1,0D0/BETA 


HlEy 

1C 

HLEW 

20 

MLEW 

30 

MLEW 

40 

MLEW 

50 

HLEW 

6C 

HLEW 

70 

MLEW 

80 

MLEW 

00 

MLEW 

IOC 

MLEW 

110 

MLEW 

120 

MLEW 

130 

MLEW 

14C 

MLEW 

15C 

MLEW 

160 

MLEW 

17C 

MLEW 

18C 

MLEW 

19C 

MLEW 

20C 

MLEW 

21C 

MLEW 

220 

MLEW 

230 

MLEW 

24C 

MLEW 

250 

MLEW 

260 

MLEW 

270 

MLEW 

280 

MLEW 

29C 

MLEW 

3  0C 

MLEW 

31C 

MLEW 

32C 

MLEW 

3  3C 

HLEW 

340 

MLEW 

350 

MLEW 

36C 

MLEW 

370 

MLEW 

36C 

MLEW 

390 

MLEW 

400 

MLEW 

410 

MLEW 

420 

MLEW 

430 

MLEW 

440 

MLEW 

450 

MLEW 

460 

MLEW 

470 

MLEW 

480 

MLEW 

490 

MLEW 

500 

MLEW 

510 

MLEW 

52C 

MLEW 

53C 

MLEW 

540 

MLEW 

5  50 

MLEW 

560 

MLEW 

570 
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OIFF-A-D  «LEW  580 

LAMBDA-TOTFAL/DENOf'  ”LEW  5<50 

WRITE  (6»12C)  BETA,D»01Ff»LAHBDA  MLEW  6CC 

ABDI«DABS(OIFF)  MLEW  61C 

If  (ABDI.LE.tPSlLN)  GO  TO  5C  MLEW  620 

IF  (DIFF.GT.O.GDO)  GO  TO  AO  MLEW  63C 

NPDIFF-J  MLEW  6AC 

BETA-BETA-APJ  MLEW  650 

If  (RETA.LE.l.OC-15)  GO  TO  60  MLEW  66C 

60  TO  10  MLEW  670 

C - MLEW  680 

C  BETA  TOO  SMALL  --  DECREASE  STEP  SIZE  AND  OSE  PREVIOUS  BETA,  MLEW  60C 

C - MLEW  7v.O 

AO  PETA»ReTA+ADJ  MLEW  71C 

NPDIFF*NPDIFF+1  MLEW  720 

IF  (KPDIFF.GT.IC)  GO  TO  70  MLEW  730 

ADJ-l.OO-l+AOJ  MLEW  7AC 

IF  (DABS(DL0G10{CABS(8ETA))-0L0G10(ACJ)).GT,1.50+1>  GO  TO  80  mlew  75^ 

BETA-BbTA-ADJ  MLEW  760 

GO  TO  10  mlew  77C 

C  ~  - - - - ________ - - - - ___ - ___  mlew  780 

C  EPSILON  TOLERANCE  MET.  MLEw  7G0 

t - mlew  800 

5C  WRITE  <6»1A0)  BE T A, L AMBD A» ABD I» E P S IL N# AO J  MLEW  810 

GO  TO  IJD  mlew  82C 

C  MLEW  830 

60  IZEBO-IZERO+i  MLEW  8A0 

IF  (IZEPO.LE.e)  GO  TO  AO  MLEW  850 

WRITE  (6>160)  MLEW  86C 

GO  TO  90  MLEW  870 

70  WRITE  (6>i7o)  MLEW  86C 

GO  TO  9o  ■  ■  mlew  890 

8D  WRITE  (6/180)  mlEw  9C0 

90  NOGDOO«.TRUE.  MLEW  910 

C  mlew  920 

100  RETURN  MLEW  933 

C  mlew  9Aj 

C  mlew  95C 

110  FORMAT  (65H1  ESTIMATION  OF  THE  PARAMETERS  OF  THE  WEIBULL  INTENS ITYMLEW  96C 

1  FUNCTION./)  mlew  970 

120  FORMAT  ( IH  > A ( 1PD16. 9, AX ) )  MLEW  96C 

13C  FORMAT  ( 3 5HCC CNS T ANT  NOT  INVOLVING  BETAt  A  «  /1PD15,7/)  MLEW  99t 

lAG  format  (30H1THE  FINAL  ESTIMATE  OF  BETA  ■  / 1PD15. 7//32H  THE  FINAL  EMLEwlOCO 
ISTIMaTE  of  LAMBDA  »  / 1  MO  1 5 . 7 / /16H  CONVERGENCE  TO  ,1PD15.7/30H  WHICMLEwlOlO 
2H  IS  LESS  THAN  EPSILON  «  » IP D 1 5 , 7 / /2 AH  THE  FINAL  STEP  SIZE  IS  /1PDMLEW1C23 
315.7/)  MLEW1C3Q 

150  format  {lH0/i6H0  ESTIMATED  eETA/5X,15H  FUNCTION  0 ( B • ) » 8 X, lOH  A  -  DMLEWlOAC 

1(B')»6X/17H  ESTIMATED  LAMBDA/)  MLEWlCiC 

160  FORMAT  (AlHO  MLEW  ERROR  —  BETA  LESS  THAN  0. 00000001)  MLEW1C6C 

176  FORMAT  (53Hri  MLEW  ERROR  —  INITIAL  ESTIMATE  Of  BETA  IS  TOO  SMALL)  MLEW1C7C 

180  FORMAT  (50Hf  MLEW  ERROR  ~  STEP-SIZE  HAS  BECOME  INS IGNIF I CANT/2CH  MLEwluBO 

IBETA  NOT  RESOLVABLE/)  MLEW1090 

END  MLEWllOO 
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SUBROUTINE  PRNTX  (  T-FAIL*  IB*  IE  I 

COHHENT  THIS  SUBROUTINE  PRINTS  THE  FAILURES  •XUB)*  THROUGH 

C  'XtlEI*. 

C 

DOUBLE  PRECISION  X(NFAIL) 

DOUBLE  PRECISION  XTEPP<5». 

DATA  L  /5/ 

C 

IT«C 

00  20  I«1B,IE 
1T«1T+1 

XTEMP( IT»»OABS<X(I)  I 
IF  (l.EO.IEI  GO  TO  10 

IF  (II.LT.LI  GO  TO 

Ij  WRITE  (6*30)  (XTEMP(J)*J»1» IT) 

IT*C 

20  CONTINUE 
C 

RE  TURN 
C 

3C  FORTAT  ( IH  * 30X *  X P5D10 . 3 ) 

END 


PRNTX  10 
PRNTX  2C 
PRNTX  30 
PRNTX  AO 
PRNTX  50 
PRNTX  60 
PRNTX  70 
PRNTX  80 
PRNTX  90 
PRNTXIOC 
PRNTXllO 
PPNTX120 
PRNTX130 
PRNTXIAO 
PRNTX150 
PPNTXI63 
PRNTX170 
PRNTX180 
PRNTX190 
PPNTX200 
PRNTX210 
PRNTX220 
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SUBPaUTINF  REPEAT  ( X»C » Nf AI L» NSYS* K )  REPET  10 

COMMENT  THIS  SUBROUTINE  REPEATS  THE  INPUT  DATA,  »EPET  20 

(■  REPET  30 

aCUBLE  PRECISION  X(NEAIL)»  C(NSYS»2)  REPET  AO 

C  REPET  50 

WRITE  (6,20)  REPET  60 

ip.[  REPET  70 

C  REPET  80 

00  10  J«1,K  REPET  90 

WRITE  (6,30)  C(  J,1),C( J»2)  REPETlOO 

IB.IP*1  PEPETllO 

CALI  LOCIP  ( X,NFAIL,M, IP)  REPET120 

Ij-.jp  REPET130 

CALL  PRNTX  ( X,Nf AIL, IB, IE  )  REPETIAO 

lu  CONTINUE  REPET150 

Q  RbPETi60 

RETURN  ■  RFPE117C 

C  REPET18C 

2C  FORMAT  < IHO / 3 A  HO S Y S T t «  STARTING  AND  ENDING  T I M E S . / IHO , 3 5X , lOH  FAILREPET190 
luRES./)  REPET200 

3C  FORMAT  (IH  ,1P2D12.3)  REPET210 

CKO  •  REPET22C 
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SUBROUTINE  SOBTX  <X»N» 

COMMENT  THIS  SUBROUTINE  SORTS  THt 

C 

DOUBLE  PRECISION  X(NJ 
M.N 

10  M-M/? 

IF  (M.EQ.O)  RETURN 
K»N-M+1 
J  ■  1 
20  I«J 
3C  L«I+M 

IF  (X<  n.GT.XIL)  )  GO  TO  5C 
40  J*J+1 

IF  (J-K»  20»10»10 
bO  T»X(L) 

x(L)*x(n 

x(n-T 

I»I-M 

IF  (1)  40/40»30 

ENO 


SORTX  10 

VECTOR  X  INTO  INCREASING  ORDER.  SORTX  20 

SORTX  30 
SORTX  40 
SORTX  50 
SORTX  60 
SORTX  70 
SORTX  80 
SORTX  9C 
SORTXIOC 
SORTXllC 
S0RTX120 
S0RTX13C 
S0RTX140 
SDRTX150 
SORTXIBO 
S0RTX17O 
S0RTX18C 
S0RTX19C 
SORTX200 
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SUBROUTINE  TRANSF  C Xf TO# IB# IBT# IE T#NF AI L ) 

COMMENT  THIS  SUBROUTINE  TRANSFORMS  THE  FAILURE  TIMES, 

C 

DOUBLE  PRECISION  X(NfAIL)#  TO 

C 

IF  CIcT.LE.O)  RETURN 
J»1B 

DO  10  I«IBT#IET 
J)/TQ 

J*J  +  i 

10  CONTINUE 
RETURN 
END 


TRANF  1C 
TRANF  20 
TRANF  30 
TRANF  40 
TRANF  50 
TRANF  60 
TRANF  70 
TRANF  eC 
TRANF  40 
TRANFIOO 
TRANFllO 
TRANF120 
TRANF13C 
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SU8R0UIINE  INREAL  (CHANNL»OESTIN» 

CQhhENT  FREE  FIELD  READ* 

C  lA  FORTRAN  TRANSLATION  OF  ACM  ALGORITHM  239. » 

r  EACH  CALL  OF  THIS  SUBROUTINE  MILL  READ  ONE  REAL  NUMBER 

C  FROM  UNIT  'CHANNL*#  CONVERT  IT»  AND  STORE  IT  IN  'OESTIN*. 


INTEGER  CHANNL 
DOUBLE  PRECISION  DESTIN 
PEAL  SIG»  FP»  D 
INTEGER  ESI6»  EP»  IP»  CH 
INTEGER  CHAR*  UNSGND 


S IG-l.O 

EP-0 

FP«C 


1C  CH«CHAR(CHANNL) 

SUPPRESS  INITIAL  BLANKS. 
If  (CH.E0.14J  GO  TO  10 

12  >  and  11  •  *-'• 

IE  (CH.NE.12)  GO  TO  20 
CH-CHAR (CHANNL ) 

GO  TO  30 

2C  IF  (CH.NE.ll)  GO  TO  30 
SIG*"l.v 
CH»CHAR(CHANNL) 

30  CONTINUE 

IF  (CH.GT.l'f)  GO  TO  70 
IF  (CH.Ct.lOI  GO  TO  <iO 
IP«tNSGND(CHANNL*CH> 

GO  TO  50 
40  CONTINUE 
I  P«C 

50  CONTINUE 

IF  (CH.NE.13)  GO  TO  130 

CH«CHAR(CHANNL  J 

FP«0 

IF  (CH.GE.IC)  GO  TO  100 

0*0.1 

60  FP*FP+f LOAT(CH)*D 
0*0*G. I 

CH*CHAR (CHANNL) 

IF  (CH.LT.IO)  GO  TO  60 
GC  TO  IfO 
7'  CONTINUE 

IF  (CH.NE.13)  GO  TO  80 
IP»1 

GO  TO  90 
60  CONTINUE 

WRITE  (6*180) 

STOP 

90  CONTINUE 
ICO  CONTINUE 

IF  (CH.NE.13)  GO  TO  16G 


INREL  10 

INREL  20 

INREL  30 

INREL  40 

INREL  50 

INREL  60 

INREL  70 

INREL  80 

INREL  9C 

INRELIOC 

INRELllO 

INREL120 

INREL130 

INREL140 

INREL150 

INREL160 

INREL17C 

INREL180 

INREL190 

INREL200 

INREL21t' 

INREL220 

INREL230 

INREL240 

INREL250 

INReL260 

INReL270 

INREL2eO 

INREL29C 

INREL30C 

INREL310 

INREL32C 

INREL330 

INREL343 

INREL350 

INPEL360 

INREL370 

INREL380 

INREL390 

INREL400 

INREL410 

INREL420 

INR6L430 

INREL440 

INREL450 

INREL460 

INREU470 

INREL4&0 

INREL490 

INREL509 

INREL51C 

INREL520 

INREL530 

INREL540 

INREL55C* 

INREL560 

INREL57C 
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CH«CHAR(CHANNL) 

ESIG-1 

IF  (CH.NE.12.AND.CH.NE.14I  GO  TO  110 
CH«CHAR(CHANNL) 

GO  TO  130 
110  CONTINUE 

IF  (CH.NE.ll)  GO  TO  120 
C  NEGATIVE  EXPONENT. 

ESIG  —  1 

CH«CHAR (CHANNLl 
120  CONTINUE 
13C  CONTINUE 

IF  (CH.C-E.IO)  GO  TO  140 
£P«UNSGND(CHANNL»CH»+ESIG 
GO  TO  150 
140  CONTINUE 

WRITE  <6»190) 

STOP 

15C  CONTINUE 
16w  CONTINUE 
C 

IF  (Cri.NE.14)  GO  TO  170 
C 

0ESTIN«06L£(SIG*(FL0AT(IP)*FP)*(10.0**EP)) 

c 

et  TURN 

l7C  WRITE  (6,200) 

STOP 

C 

iaO  FORMAT  OeHC  INREAL  ERROR  —  CH  OUT  OF  RANGE.  /) 

1<?0  FORNAT  (AlHG  INREAL  ERROR  —  EXPONENT  NOT  DIGIT. 

200  FOPNaT  {53H0  INREAL  ERROR  —  NO  BLANK  FOUND  BETWEEN 

ENO 


/) 

DATA  VALUES 


INRELSeC 
INREL590 
INRE160C 
INREL610 
INPEL620 
INREL630 
INREL6A0 
INREL650 
INREL660 
INREL67C 
INREL68C 
INREL690 
INREL7C0 
INREL71Q 
INREL720 
INREL730 
INREL7A0 
INREL750 
INREL760 
INREL770 
INREL7aC 
INPEL790 
INREL800 
INRELtilO 
INREL820 
INREL830 
INREL04O 
INREL850 
INREL860 
INREL870 
INREL880 
/ ) INREL89C 
INREL9aO 
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o  o  o 


INTEGER  FUNCTION  CHAR  ICHANNL)  CHAR  10 

OMHENT  'CHAP*  RETURNS  AN  INTEGER  VALUE  FOR  THE  NEXT  CHARACTER  CHAR  20 

ON  UNIT  ‘CHANNl'. 

CHAR  40 


C 


c 


INTEGER  C 
INTEGER  CHANNL 
INTEGER  STRING(15) 

DATA  STRING(l)#  STRING(2)#  STRlNG(3)f  STRlNG<4)f 
1(6),  STRING(7),  STRING(8),  STRING(9),  STRING(IO) 
2NG(12),  STRING(13),  STRING(14),  STRING(15)  /IHO, 
34,  lrt5,  1H6,  1H7,  IHB,  lH9f  IH.,  1H-#  IH^,  1H6, 
DATA  LSTR  /15/ 

CALL  INSYML  (CHANNL, STRING, LSTR, C  ) 

IS  CHARACTER  LEGAL? 

IF  (C.LE.O)  GO  TO  10 

CHAR-C-1 

RETURN 

10  WRITE  (6,20) 

STOP 


CHAR 

CHAR 

CHAR 

STKING(5)f  STRINGCHAR 
,  STRINGdDf  STRICHAR 
iHlf  iHZf  1H3,  IHCHAR 
IH  /  CHAR 

CHAR 

CHAR 

CHAR 

CHAR 

CHAR 

CHAR 

CHAR 

CHAR 

CHAR 


20  FORHAT  (58H3  CHAR  ERROR  —  ILLEGAL  INPUT  CHARACTER.  PROGRAM  ABORTCHAR 
lED./) 

END  CHAR 


50 

60 

70 

80 

90 

100 

110 

120 

130 

140 

150 

169 

170 
180 
190 
200 
210 
22C 
230 
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SUBROUTINE  INSYML  ( lUN IT# STRING, L SIR t 1) 

COMMENT  THIS  SUBROUTINE  EXAMINES  THE  NEXT  CHARACTER  ON  ‘lUNIT* 

C  AND  DETERMINES  ITS  POSITION  NUMBER  WITHIN  THE  •STRING*  OF 

C  LENGTH  "LSTR*. 

C 

INTEGER  STRING(LSTR) 

INTEGER  BUFFER<160) 

DATA  IP  /O/ 

DATA  LENGTH  /80/ 

C 

IF  (IP.NE.O)  GO  TO  1C 
C  FILL  INPUT  BUFFER. 

READ  (IUNIT,20)  < BUFF ER ( J > # J«l# L ENGTH ) 

1C  CONTINUE 
IP*IP+1 
IC«IP 

C  MATCH  THE  CHARACTER. 

CALL  MATCH  ( BUFFER < IC ), STR I NG, L STR , 1) 

C  IF  POINTER  'IP*  HAS  REACHED  THE  END  OF  A  LINE#  RESET  IT. 

IF  (IP. EQ. LENGTH)  IP«0 

C 

RETURN 

C 

20  FORMAT  (128A1) 

END 


INSYL  10 
INSYL  20 
INSYL  30 
INSYL  ^0 
INSYL  50 
INSYL  60 
INSYL  70 
INSYL  80 
INSYL  90 
INSYLIOO 
INSYLIIO 
INSYL120 
IN5YL130 
INSYL140 
INSYL150 
INSYL160 
INSYL170 
INSYL180 
INSYL190 
INSYL20^' 
INSYL210 
INSYL220 
INSYL230 
INSYL240 
INSYL25^’ 
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SU8K0UTINE  MATCH  ( CHAR# STR ING# LSTR^ I P ) 

COMMENT  THIS  SUBROUTINE  FINDS  THE  POSITION  'IP*  OF  •CHAR*  IN 

C  'STRING*  WHICH  HAS  A  LENGTH  OF  'LSTR*. 

C 

INTEGER  CHAR,  STRING(LSTR> 

IP*r 

00  10  1«1,LSTR 

IF  ($TRING(n.N£,CHAR>  GO  TO  1C 
IP«I 

GO  TO  2C 
ID  CONTINUE 

WRITE  <6#30)  CHAR 
20  RETURN 
C 


MATCH  10 
MATCH  20 
MATCH  30 
MATCH  40 
MATCH  50 
MATCH  60 
MATCH  70 
MATCH  8C 
MATCH  90 
MATCHIOO 
MATCHllO 
MATCH120 
MATCH130 
MATCH14C 


30  FORMAT  (31H0  MATCH  ERROR  —  THE  CHARACTER  ,A1,16H  IS  NOT  MATCHED. /MATCH15C 


1) 


MATCHlbO 


END 


MATCH17C 
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INTEGEB  FUNCTION  UNS6ND  (CH4NNl,CH»  UNSG  10 

COMMENT  THIS  FUNCTION  RETURNS  THE  NEXT  UNSIGNED  INTEGER  FROM  UNSG  20 

C  'CHANNL'.  UNSG  30 

C  UNSG 

INTEGER  CHANNL*  CH  UNSG  50 

INTEGER  CHAR  UNSG  60 

INTEGER  U  UNSG  70 

g.O  UNSG  eo 

10  U»1C*U+CH  UNSG  90 

CH«CHAR(CHANNL)  UNSG  100 

IF  (CH.LT.IO)  GO  TO  10  UNSG  110 

UNSGNO*U  UNSG  120 

return  UNSG  130 

EE)0  UNSG  lAO 


Next  page  is  blank. 
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APPENDIX  B 
INPUT  FOR  TEST  CASES 


Next  page  is  blank. 
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INPUT  DATA  FCP  TEST  CASE 


o 


CO 


O 

O' 

• 

• 

CO 

CO 

CO 

o 

• 

C- 

rH 

r-4 

r^ 

r- 

• 

• 

O' 

O' 

o 

rH 

00 

o 

• 

• 

o 

O' 

rO 

rH 

Cr 

Co 

• 

• 

rH 

iH 

• 

ro 

O 

t-H 

O 

O' 

• 

c* 

• 

CJ 

O' 

1^ 

in 

rH 

• 

• 

OCl 

r-i 

Xf 

rH 

OD 

CO 

O 

OD 

rH 

• 

• 

rH 

CO 

ro 

f'- 

rH 

C- 

fw 

>r 

O' 

• 

• 

O' 

C- 

fNJ 

rH 

• 

o 

rH 

• 

• 

CO 

sj- 

>o 

• 

• 

00 

CO 

T-H 

of 

u-^ 

O' 

rH 

u 

O' 

r- 

• 

r~4 

ro 

O' 

O' 

o 

O' 

r-  1 

HO 

O' 

• 

• 

1 

or 

« 

1 

00 

r- 

iTS 

• 

• 

ro 

• 

00 

O' 

• 

OD  O' 

vO 

C- 

• 

cr 

O' 

<) 

CO 

CD 

1  O' 

CO 

CO 

CO 

• 

a 

O' 

rH 

CO 

1 

O' 

O' 

1 

O' 

vi) 

O' 

rH 

oo 

>D 

X 

CO  X 

• 

1 

O' 

U'> 

• 

f'-  O 

• 

•  • 

‘O 

O' 

• 

00 

• 

OC-i 

CO 

or  o 

u.\ 

• 

O' 

1 

o* 

04 

a- 

O' 

• 

a> 

u:^ 

• 

X 

rH 

CO 

1 

a\ 

Co 

O') 

4J 

O' 

O' 

rH 

CO 

JD 

U'i 

• 

O' 

• 

vO 

'O- 

m 

m 

• 

•  O 

CD 

O' 

• 

• 

CO 

^  >r 

ro 

O'  • 

• 

1 

• 

ur\ 

r- 

OD 

O' 

CO 

rH  O 

c-- 

jr\ 

•j- 

r-t 

U". 

• 

a;i 

r'- 

rH 

IL^ 

X-' 

(O 

rH 

C- 

or 

• 

c- 

• 

CO  U'  T 

CO 

>D 

• 

O' 

a.> 

rH 

CO 

'D  'T 

• 

• 

O' 

m 

• 

• 

m 

rH 

r- 

o 

'JD  CO 

<r 

• 

sO 

ro 

• 

o 

rH 

rH  • 

r- 

• 

ijr> 

r*^ 

r-» 

ro 

CO 

•  r-H 

rH 

vO 

O' 

rH 

U'N  • 

sO 

tn 

<r 

00 

Ln 

• 

^  vO 

'0- 

sD  rH 

• 

• 

r- 

sO 

O' 

ITi 

• 

« 

CO 

• 

• 

• 

ro 

• 

O' 

m 

CO 

• 

CO 

nT 

CO 

^x> 

o> 

rH 

oo 

o 

rH 

OD 

OD 

rH 

• 

• 

rH 

CO 

Co 

• 

• 

in 

fVl 

o 

O' 

ro 

rH  O 

>r  m 

>r 

O 

• 

CD 

O' 

rH 

u> 

• 

•  fvi 

•  O' 

• 

vN4 

O 

• 

CO 

o  >r 

O' 

O  rH 

vD  rH 

in 

ro 

• 

• 

• 

• 

sO 

o 

%r 

ro 

o 

tfT 

ro 

O 

CO 

CO 

• 

CO 

CD 

CD 

sT 

O 

• 

• 

• 

o 

CO 

• 

• 

• 

• 

o 

• 

r-4 

00 

rH 

o 

• 

O  O' 

o 

m 

• 

O 

CO 

O' 

00 

O' 

•  CO 

O' 

O' 

O' 

O 

'4' 

o 

• 

CO  r-l 

>r 

rH 

O' 

O' 

rH 

rH 

rH 

• 

fO 

rH 

1 

r- 

o 

in 

-r 

• 

CO 

O' 

O' 

O' 

o 

in  csj 

a  CO 

• 

• 

o 

o 

O'  o 

O' 

O' 

CO 

U 

u  •  • 

CO  • 

HO  • 

O  HO 

•  • 

• 

O'  • 

O' 

O' 

O' 

• 

cj  •  ^  m 

O  O' 

m  O' 

•  m 

O  O' 
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tILEW  ERROR  —  BETA  LESS  THAN  O.UOCOOOOl 
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CONSTANT  NTT  INVOLVING  BETA:  A  ■  3.6731896D+C0 


<L 

fVJ 

O 

OO 

(O 

CTI 

o 

4VJ 

r-4 

•-4 

r4 

r-4 

fH 

t-4 

r-4 

t-^ 

t-4 

t-4 

t-4 

r-4 

t-4 

t-H 

t-4 

t-4 

r-4 

t-4 

r-4 

1-4 

o 

OJ 

<M 

rH 

r-4 

iH 

1-4 

o 

o 

o 

Q 

o 

O 

O 

o 

o 

o 

o 

o 

o 

O 

o 

o 

O 

o 

o 

o 

O 

o 

o 

O 

o 

CD 

I 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

• 

1 

1 

1 

I 

1 

1 

1 

1 

1 

1 

1 

S. 

O 

O 

o 

O 

O 

o 

o 

o 

o 

o 

o 

o 

O 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

< 

to 

40 

O' 

fH 

lO 

-O' 

ro 

rsi 

40 

o 

to 

yO 

to 

>o 

ro 

r-4 

40 

CO 

yO 

yf 

O 

r- 

CO 

m 

o 

r-l 

o 

OJ 

rn 

fSJ 

CO 

O" 

O' 

40 

r- 

r- 

csj 

r- 

O 

•o 

m 

X 

r- 

fO 

O' 

O 

yT 

4- 

o- 

h- 

yO 

O' 

CO 

<v 

o 

(NJ 

o 

O 

K 

o 

Ctl 

o 

QO 

r-4 

fvj 

<NJ 

X 

ct> 

f'- 

r*- 

<NJ 

o 

CO 

O 

>D 

o 

<r 

to 

CO 

O 

Ul 

t-4 

o 

«Si 

•H 

r- 

m 

r*- 

1-4 

m 

oo 

vT 

lO 

40 

4i> 

40 

to 

IM 

r-4 

r- 

o 

o 

t-4 

CO 

o 

-O' 

CO 

o 

o 

yt 

o 

ao 

t-4 

lO 

LU 

o 

m 

iH 

m 

o 

cn 

O' 

h- 

O 

r-l 

O 

1-4 

40 

>o 

O 

(VI 

4* 

CO 

fVi 

c- 

CO 

o 

<NJ 

o 

o 

CO 

r- 

yO 

c 

CO 

o 

or 

40 

X' 

m 

f-4 

>0 

U'^ 

t-4 

O 

O 

CO 

fO 

'C 

00 

-r 

o 

CO 

ro 

40 

oo 

<M 

t-4 

'T 

CO 

o 

on 

< 

m 

O 

rh 

1-4 

r* 

f'- 

on 

>o 

o 

X 

O 

O 

tF^ 

t-4 

1-4 

ro 

to. 

o 

OO 

r-l 

>r 

o- 

r-l 

O 

CO 

O 

o 

O 

i:' 

(M 

o 

CO 

O 

O 

uv 

fVJ 

o 

40 

CO 

1-4 

o 

O 

t-4 

m 

o- 

lO 

X 

o 

CO 

40 

CO 

CO 

CO 

fO 

CO 

40 

O' 

CO 

iH 

«M 

<NJ 

m 

>o 

•H 

»-4 

o 

r- 

o 

o 

r-4 

CO 

44\ 

o 

r-4 

CO 

UA 

CO 

CO 

■4- 

>r 

yf 

'O- 

>r 

'T 

t/> 

>— 4 

INi 

40 

fVI 

cu 

r-4 

m 

r-4 

fM 

to 

m 

CO 

m 

CO 

fn 

yT 

yt 

-r 

nT 

-r 

•4' 

>r 

^r 

yf 

o 

ID 


o 

o 

o 

o 

O 

o 

o 

t-4 

1-4 

t-4 

f— 1 

CO 

-4 

fH 

tH 

CO 

CJ 

Cd 

CJ 

CO 

m 

CO 

ro 

CO 

to 

•r 

ifi 

X 

u 

o 

V' 

o 

«.> 

o 

<-> 

CP 

o 

O 

o 

:l> 

o 

o 

o 

o 

o 

o 

O 

o 

o 

'D 

'D 

•4- 

+ 

♦ 

+ 

-f 

■t* 

4- 

•4- 

1 

1 

1 

1 

1 

1 

1 

1 

1 

I 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

«-k 

o 

o 

Q 

C. 

o 

o 

o 

O 

O 

O 

o 

o 

O 

o 

Ct 

o 

o 

o 

o 

o 

o 

o 

Q 

o 

o 

o 

o 

o 

o 

Cj 

O 

O 

m 

40 

00 

CO 

CO 

CO 

CO 

40 

yf 

rH 

o 

X 

O' 

•L> 

o 

O' 

X 

y^ 

r- 

CO 

o 

N- 

O' 

to 

tH 

c- 

X 

o 

CO 

o 

O' 

ro 

fO 

CJ 

>r 

X 

-r 

40 

O* 

X 

o 

to 

X 

t-4 

lO 

X 

•r 

CO 

fH 

rH 

to 

CO 

ro 

fO 

m 

X 

lO 

O 

•O’ 

*rt 

r- 

to 

yO 

CO 

c- 

CO 

CO 

X 

CO 

OO 

•O' 

to 

CO 

1-4 

C- 

tH 

O' 

o* 

CO 

•O’ 

to 

O' 

C4 

U 

X 

tH 

O 

ro 

o 

X 

o 

'4' 

ry 

yD 

CO 

co 

o 

sC 

X 

'T 

'J' 

fO 

40 

c- 

CO 

CO 

CO 

c- 

yt 

X 

o 

cn 

tH 

CO 

O 

X 

X 

X 

o 

•0- 

O' 

o 

r- 

X 

O 

yu 

X 

CO 

-O' 

r-l 

to 

X 

4J 

40 

o 

O' 

tr 

ro 

o 

o 

4'« 

o 

X 

o 

u-t 

o 

•CD 

c- 

fv. 

CO 

1 

lO 

CJ 

Cd 

rH 

1-4 

X 

>r 

O' 

o 

trt 

CO 

1-4 

CO 

o 

1-4 

o 

40 

o- 

to 

iH 

CO 

00 

O'l 

X 

o 

•r 

X 

X 

O' 

O' 

tH 

O' 

ao 

X 

Ol 

X 

X 

QO 

X 

r- 

X 

CO 

O' 

tO 

CO 

CO 

CO 

CO 

X 

X 

CO 

to 

m 

r- 

X 

X 

CO 

O' 

O' 

< 

O' 

X 

0- 

40 

CO 

O' 

yC 

O' 

c- 

sO 

c- 

O' 

VC 

X 

to 

yO 

CO 

CO 

tH 

to 

CO 

X 

r- 

rH 

-r 

CJ 

O' 

X 

X 

'T 

yf 

>r 

CO 

m 

o 

a 

CO 

•r 

X 

X 

CO 

o 

00 

O' 

to 

o 

in 

O' 

■•» 

X 

co 

CO 

CO 

c- 

X 

fH 

X 

— 1 

*— 4 

t—4 

rH 

tr. 

yS 

ro 

rH 

>o 

tH 

fH 

X 

to 

CO 

tH 

c- 

If  : 

CO 

rO 

CO 

a 

iH 

in 

X 

u^ 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

O 

O 

o 

o 

o 

o 

o 

c- 

c> 

o 

o 

O 

o 

O 

o 

o 

o 

o 

o 

o 

rv 

o 

o 

o 

o 

o 

o 

r  ^ 

o 

o 

o 

o 

rj 

o 

o 

o 

o 

Q 

o 

o 

o 

O 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

b 

c> 

m 

+ 

+ 

•¥ 

+ 

•I- 

•f 

•f 

+ 

•f 

4' 

•f 

♦ 

+ 

•f 

•f 

4- 

+ 

+ 

+ 

+ 

+ 

-4- 

-f 

4- 

•♦• 

+ 

+ 

CD 

o 

Cl 

o 

O 

o 

o 

c- 

o 

o 

c- 

o 

C 

o 

o 

o 

o 

c 

o 

o 

O 

o 

o 

o 

Q 

O 

o 

o 

o 

o 

o 

o 

Cj 

w 

nT 

X 

CM 

CO 

fH 

•O' 

•H 

X 

X 

X 

c- 

O' 

iTi 

X 

c- 

o 

o 

X 

X 

CO 

l/l 

o 

CM 

X 

o 

o 

p- 

CM 

X 

X 

cc 

Cj 

m 

•O' 

X 

CO 

X 

O' 

tH 

O' 

X 

X 

O' 

O' 

f"- 

X 

O' 

CO 

CO 

O' 

X 

to 

tH 

CO 

to 

o 

X 

o 

X 

p- 

•J- 

fH 

X 

CO 

CM 

c 

X 

O' 

O' 

ac 

•f- 

C' 

o 

rH 

to 

X 

O' 

O' 

c- 

X 

CJ 

X 

o 

O' 

tH 

X 

X 

X 

rH 

r- 

*o 

o 

tr 

S' 

sT 

o 

X 

ro 

X 

X 

u't 

nT 

X 

tn 

X 

•r 

O' 

X 

X 

CJ 

X 

>r 

•J- 

o 

X 

to 

U'\ 

r- 

tH 

CM 

X 

X 

rH 

U*'* 

o 

X 

c- 

O' 

\o^ 

m 

CM 

X 

o 

•J’ 

-O' 

tH 

»H 

O' 

o 

O' 

•r 

CM 

rH 

n- 

r- 

X 

•r 

r4 

c- 

CM 

c- 

r- 

iH 

or 

CM 

a 

t— t 

c- 

yT 

X 

X 

X 

X 

X 

tH 

CM 

<;> 

O' 

o 

vT 

CM 

h- 

X 

yf 

X 

in 

•0- 

X' 

o 

tH 

X 

X 

o 

CM 

o 

CO 

CM 

tH 

4- 

CO 

rH 

c- 

O' 

to 

tH 

tH 

X 

tH 

rH 

O' 

fH 

CM 

•^^ 

tH 

o 

X 

to 

CO 

X 

X 

X 

rH 

X 

X 

X 

o 

X 

X 

X 

X 

X 

tD 

X 

nT 

CO 

o 

r- 

X 

X 

c- 

U'4 

CO 

nT 

n 

CM 

o 

U  1 

X 

o 

X 

to 

X 

r- 

c- 

p- 

p- 

p^ 

p- 

r- 

c- 

2. 

r4 

fH 

tH 

rH 

rH 

o 

c,> 

O' 

C- 

CM 

tH 

o 

X 

X 

X 

X 

c- 

c- 

r- 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

u. 

lO 

u'A 

U> 

to 

in 

to 

•O’ 

>1“ 

•O’ 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

rH 

i:> 

o 

o 

o 

c* 

o 

o 

o 

o 

r4 

tH 

-4 

rH 

t-4 

r-4 

r4 

rH 

rH 

rH 

rH 

iH 

tH 

»H 

tH 

tH 

iH 

fH 

iH 

r-4 

1^ 

<3, 

o 

o 

o 

o 

o 

c 

o 

o 

o 

o 

«*rC 

c 

o 

o 

o 

o 

o 

o 

o 

o 

c 

o 

o 

z> 

o 

o 

o 

o 

o 

r> 

o 

K- 

•p 

•f 

+ 

•p 

•p 

•p 

•p 

4^ 

■p 

■p 

1 

1 

1 

1 

1 

1 

1 

1 

• 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

I 

1 

1 

X 

o 

c. 

o 

o 

o 

Cl 

o 

o 

c 

C'. 

o 

a 

o 

c 

Cl 

c 

c 

c 

c 

o 

o 

o 

Q 

c 

Cj 

c 

Q 

c 

o 

o 

c 

CD 

o 

o 

o 

o 

r> 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

D 

o 

o 

o 

c> 

o 

X 

o 

Q 

o 

o 

o 

o 

? 

u 

o 

o 

o 

o 

o 

tD 

o 

o 

o 

o 

o 

6 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

w 

ID 

o 

o 

CD 

o 

o 

o 

o 

o 

o 

ID 

(D 

o 

o 

o 

o 

o 

o 

CD 

o 

o 

o 

<D 

o 

o 

Q 

o 

CD 

LU 

ID 

o 

o 

o 

o 

c> 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

1— 

Q 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

Q 

o 

o 

< 

o 

o 

o 

CD 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

O' 

X 

c- 

z: 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

c 

o 

o 

o 

o 

o 

o 

o 

o 

o 

O' 

cr: 

X 

X 

X 

♦—4 

o 

o 

o 

o 

o 

CD 

a 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

O' 

X 

f"- 

X 

to 

in 

U  1 

U  1 

in 

to 

X 

UJ 

o 

o 

o 

o 

o 

CD 

o 

o 

o 

o 

o 

o 

& 

o 

O' 

X 

X 

to 

•r 

X 

CJ 

CM 

CJ 

CO 

CJ 

CO 

CO 

CO 

CO 

CM 

CM 

rH 

O' 

X  pw 

X 

to 

X 

CO  IH  O' 

X 

X 

X 

X 

XXX 

X 

X 

X 

X 

X 

•c 

X 

X 

X 

X 

X 

X 

X 

62 


THE  FINAL  --.STIMATF  OF  BETA  =  6 . 2  E  87C0DD-O  i 


o 


1 

in 

Q 

o 

1 

Q 

rn 

O 

TO 

o 

o 

o 

in 

ro 

o 

■-3 

o 

1 

• 

c> 

O 

nT 

o 

o 

• 

o 

o 

o 

H 

sO 

r> 

o 

O 

< 

1  n 

€*1 

O 

o 

'•  • 

03 

TTj  :2 

f— i 

x: 

o  o 

< 

fM  ^ 

«j 

r— 1  i-H 

O'  tyO 

yy 

vO  a- 

D 

a’  'U 

• 

'Li 

OJ 

U'l  -^1 

hxj 

♦— 

< 

F-H 

< 

X 

tn 

r. 

♦— 

1-^ 

o 

a. 

K- 

h-  yy 

LU 

L-O 

i/y 

I-’- 

■_u 

Ui  Lit 

1/3 

O  -J 

-J 

< 

iXi  t/3 

-< 

O 

z 

QC 

F-H 

UL 

ay  X 

U. 

:>  o 

Oi 

z 

UJ 

X 

O  X 

X 

♦- 

3 

H- 

63 


TOTAL  NUMBtR  OF  SYSTEMS 


<M  f-l 

O  O 

4-  ♦ 

o  o 
o  o 
^  o 

CD  CO 

•  • 

IM  »-♦ 
O  O 
4>  4* 

Q  O 

o  o 
o 

\n 

m  • 

fM  CM 
O  O 
4-  + 

O  O 

o  o 

>r 

sO 

•  • 

1^  ir\ 


CM  r4 

rH 

O  O 

O 

4-  1 

4* 

o  o 

Q 

<.;>  o 

O 

• 

'O  o 

rH 

00 

00  CM 

O' 

UJ 

•  • 

• 

oc 

M*  rH 

O' 

r-l  CM 

rH 

t-4 

o  o 

o 

< 

4-  ♦ 

4 

u. 

O  Q 

o 

o  o 

r> 

c>  CM 

O 

LCO  O" 

O' 

•  • 

• 

'J'  M* 

CD 

u 


o 


CM 

o 

rs 

+ 

ai 

o 

c- 

O 

1C 

o 

CO 

o 

2: 


o: 

< 

t- 

t/> 


s: 

1X1 

♦— 

</) 

V 

*/> 


65 


ESTIMATED  BETA  FUNCTION  0(8»)  A  -  D(B»)  ESTIMATED  LAMBDA 


00 

m 

o 

00 

kO 

CVJ 

e^j 

fv) 

tH 

r-l 

rH 

tH 

r-l 

tH 

tH 

fH 

fH 

tH 

iH 

iH 

fH 

tH 

fH 

fH 

tH 

rH 

rH 

tH 

tH 

tH 

rH 

fH 

tH 

rH 

tH 

<VJ 

CNJ 

cvi 

M 

i-i 

o 

o 

O 

O 

o 

O 

O 

o 

o 

o 

O 

o 

o 

o 

o 

O 

o 

O 

O 

o 

o 

Q 

o 

O 

o 

O 

o 

O 

O 

o 

O 

O 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

I 

t 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

o 

Q 

o 

O 

o 

o 

O 

o 

o 

o 

O 

O 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

O 

o 

o 

o 

o 

o 

O 

o 

o 

o 

o 

d 

o 

a 

o 

o 

O 

o 

«£/ 

fn 

00 

(> 

m 

o 

o 

O 

o 

o 

yQ 

CO 

iH 

o 

in 

o 

CM 

O' 

00 

o 

4 

O' 

r- 

CO 

00 

CD 

in 

in 

o 

00 

O' 

4> 

00 

00 

4} 

m 

4 

tH 

r- 

o* 

r- 

o 

o 

o 

Q 

CT* 

fH 

O' 

oo 

o 

yO 

m 

oo 

CM 

o 

in 

»H 

00 

r'- 

CO 

tH 

o* 

o 

in 

CM 

CM 

in 

o 

00 

O' 

m 

O 

o 

r- 

fO 

CD 

CO 

lO 

o 

o 

o 

m 

00 

CO 

O' 

yO 

iH 

tn 

r- 

yO 

4- 

O 

O' 

in 

00 

4) 

CM 

CM 

O' 

4 

CM 

O' 

CM 

CO 

CM 

o 

m 

O' 

CM 

fCl 

fw 

fw 

tn 

o 

(T' 

f*- 

o 

o 

Q 

oo 

tH 

rH 

00 

O 

CO 

oo 

CO 

yO 

oo 

o 

CO 

O 

CO 

o 

in 

CM 

o 

m 

00 

tH 

in 

rH 

O 

O 

00 

O' 

m 

m 

o 

00 

tr\ 

o 

o 

>»* 

o 

>r 

CO 

'.I* 

m 

O 

in 

CM 

iH 

*H 

o 

4 

CM 

*H 

O' 

4 

oo 

f''* 

o* 

CM 

m 

O 

CM 

43 

O 

4 

<A 

ie\ 

o 

«o 

fH 

o 

Q 

00 

c> 

CO 

CO 

yQ 

yO 

CO 

fH 

00 

CM 

CM 

r- 

rO 

4) 

O 

yC 

CO 

tH 

in 

CM 

O 

fr 

4 

CM 

O' 

O'* 

4 

r- 

o 

OJ 

rsi 

(M 

a 

m 

lO 

Q 

o 

o 

>r 

>r 

CM 

yT 

in 

h- 

o 

r- 

1^ 

tH 

tH 

cu 

CM 

4 

»H 

to 

in 

00 

fH 

4 

4) 

O' 

CM 

4 

C' 

tn 

(NJ 

>o 

rj 

r* 

o 

C\J 

r-l 

tH 

yt 

CM 

OD 

in 

4* 

CO 

4- 

f'- 

rH 

O 

4 

O 

CM 

4 

o 

rr 

CO 

oc 

OO 

OO 

O' 

O 

O’ 

o 

cr 

r- 

4f\ 

cr 

CO 

00 

\n 

iH 

iH 

CM 

CO 

■J- 

aO 

CO 

o 

CM 

4 

n- 

O' 

CM 

O' 

O 

O 

Q 

o 

o 

o 

O 

o 

o 

O 

• 

o 

• 

rH 

• 

• 

r-l 

• 

CD 

• 

• 

to 

• 

• 

CO 

• 

C\J 

• 

fVJ 

• 

f-4 

• 

rH 

• 

<SJ 

in 

fH 

fVJ 

4^ 

CM 

CM 

CM 

CM 

CO 

CO 

CO 

CO 

CO 

4 

CO 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

o 

o 

O 

O 

o 

o 

o 

o 

o 

rH 

fH 

fH 

tH 

rH 

CM 

fH 

fH 

tH 

•H 

fH 

fH 

CM 

CM 

CM 

CM 

CM 

CM 

CO 

m 

4 

CO 

CO 

CO 

CO 

CO 

CO 

m 

4 

4 

o 

o 

o 

o 

o 

O 

o 

o 

CD 

CD 

CD 

O 

O 

CD 

<D 

ID 

CD 

CD 

CD 

CD 

O 

CD 

CD 

CD 

O 

O 

o 

CD 

V 

o 

OJ 

AD 

o 

'D 

iD 

CD 

o 

CD 

CD 

H* 

H 

H 

+ 

H- 

H« 

H 

H- 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

O 

o 

Q 

O 

o 

O 

Q 

Q 

O 

O 

o 

O 

Q 

o 

o 

O 

O 

O 

O 

O 

o 

Cj 

O 

O 

o 

o 

O 

O 

o 

o 

Q 

O 

o 

Q 

O 

Ci 

o 

O 

O 

o 

O' 

O 

r- 

CC> 

o 

O 

r- 

O 

GO 

p- 

00 

CO 

CM 

r^ 

43 

CO 

o 

tH 

P- 

o 

o 

4) 

tH 

K 

o 

CO 

CO 

4 

43 

4 

a 

4 

O' 

fH 

O' 

CM 

rH 

4 

CM 

o 

CM 

r- 

m 

CM 

CM 

OD 

CM 

O' 

00 

O' 

iH 

m 

tH 

in 

CD 

4 

4 

fH 

CO 

o 

P- 

tH 

fH 

in 

X 

43 

4 

P'* 

P'“ 

CM 

CM 

43 

P'“ 

m 

« 

43 

'D 

43 

n 

4> 

4 

o 

4> 

4? 

CM 

43 

iH 

o 

rH 

O' 

in 

O 

43 

tn 

in 

CO 

o 

43 

p- 

43 

4) 

o 

CM 

o 

X 

43 

4 

cr 

4 

4 

X 

CO 

X 

o 

CM 

CO 

CM 

m 

4> 

4> 

cn 

m 

o 

cn 

4> 

m 

4> 

rH 

o 

CX3 

fO 

43 

4 

m 

X 

P- 

43 

O' 

O 

to 

r- 

CO 

CM 

4 

O 

P- 

U'l 

O' 

X 

CM 

CM 

43 

4 

cr 

O' 

OO 

O' 

CO 

CM 

O' 

O' 

43 

O' 

m 

CM 

CO 

CM 

43 

CO 

(O 

m 

pr 

O' 

rH 

fH 

fH 

43 

CO 

CO 

iD 

43 

CO 

CO 

o 

X 

X 

CM 

t-H 

4 

CO 

4 

4; 

o 

O 

o 

CM 

4 

o 

CD 

O 

O' 

oc 

O' 

CM 

CM 

43 

in 

o 

»H 

4 

tH 

43 

O' 

43 

r- 

O' 

o 

o 

O' 

CO 

4 

43 

n 

CM 

X 

CM 

CM 

4 

o 

to 

4% 

r- 

CO 

o 

O 

43 

o 

o 

O' 

o 

m 

4 

o 

#H 

O' 

CO 

o 

CM 

o 

tH 

o 

O 

m 

X 

p- 

O' 

CO 

X 

P'* 

X 

iC'l 

4 

(M 

fH 

U‘i 

43 

CD 

43 

m 

CO 

tH 

00 

cn 

4 

oo 

o 

oo 

o 

rH 

CO 

o 

in 

in 

in 

4 

rH 

CO 

m 

o 

m 

O' 

4 

4 

O 

o 

4 

rH 

CM 

CO 

4 

4 

43 

P-* 

CP 

P- 

OO 

OO 

OD 

CD 

OO 

yQ 

4 

00 

43 

CO 

m 

•H 

O' 

oc 

4 

CM 

O' 

43 

CM 

CO 

O' 

O' 

4 

rH 

CO 

o 

4 

D 

4> 

CM 

X 

4 

* 

O 

• 

P*- 

• 

X 

• 

# 

• 

rH 

• 

• 

rH 

pH 

rH 

O' 

OD 

4 

X 

cH 

CM 

CM 

CM 

rH 

rH 

rH 

CP 

4 

rH 

•H 

rH 

X 

fH 

X 

X 

CM 

CM 

rH 

r-l 

X 

PJ 

1 

1 

1 

1 

1 

1 

* 

1 

1 

1 

1 

1 

1 

1 

1 

1 

* 

1 

1 

1 

1 

1 

1 

1 

I 

1 

* 

1 

f 

1 

1 

1 

1 

1 

i 

(A 

cr 

O 

r:> 

o 

o 

o 

cr* 

o 

o 

o 

o 

o 

o 

O 

O 

^D 

o 

o 

C> 

o 

O 

o 

o 

o 

D 

O 

o 

O 

D 

o 

O 

o 

o 

O 

’D 

o 

o 

o 

o 

CD 

o 

D 

Ci 

w 

o 

o 

o 

O 

CD 

o 

o 

O 

O 

o 

Q 

O 

CD 

o 

CD 

o 

o 

O 

CD 

o 

o 

vD 

o 

D 

o 

o 

;d 

o 

C'D 

O 

H- 

H- 

4> 

+ 

H- 

-f 

H- 

H- 

4- 

-f 

♦ 

♦ 

•f 

■p 

+ 

■P 

•p 

+ 

-P 

•P 

■p 

■P 

+ 

+ 

■P 

4- 

+ 

4* 

4- 

+ 

+ 

4« 

4- 

4 

4- 

♦ 

4 

O 

o 

o 

o 

O 

o 

O 

O 

D 

O 

O 

o 

o 

CD 

tD 

o 

o 

Q 

o 

o 

o 

O 

o 

O 

o 

o 

O 

O 

o 

O 

O 

CD 

O 

iD 

o 

O 

O 

C 

d 

X 

p- 

X 

iTi 

X 

X 

4 

X 

X 

P- 

X 

O' 

rH 

X 

X 

P“ 

X 

CM 

CP 

X 

4 

4 

4 

X 

4 

rH 

o 

X 

X 

4 

4 

P“ 

X 

rH 

o 

4 

4 

O 

CM 

rH 

CM 

X 

43 

CM 

CM 

cr 

CM 

CM 

CM 

CP 

43 

X 

p^ 

O' 

X 

4 

o 

p- 

P* 

o 

fH 

CM 

X 

o 

4 

4 

o 

O' 

CP 

X 

PM 

X 

4 

X 

X 

4 

P- 

43 

P- 

4 

o 

P^ 

P- 

CO 

P- 

P- 

43 

P* 

CM 

O 

p^ 

4 

p- 

4 

O' 

X 

PJ 

X 

4 

P- 

iH 

p'- 

fH 

fH 

p- 

4 

X 

p’ 

X 

X 

O' 

X 

X 

P- 

rH 

4 

a'k 

rH 

X 

rH 

»H 

O 

fH 

o 

4 

iH 

4D 

X 

4 

P- 

o 

a 

CP 

4 

P- 

tH 

4 

p- 

4 

P'^ 

X 

X 

4 

4 

X 

CM 

o 

CP 

P'~ 

4 

O 

4 

fH 

tH 

fH 

O 

tH 

4 

4 

tH 

O' 

p^ 

CM 

CP 

PM 

4 

rH 

X 

CM 

fH 

P- 

o 

PM 

PM 

tH 

X 

4 

4 

X 

p-" 

X 

X 

0 

o 

4i 

4 

43 

P- 

O' 

4> 

X 

CM 

43 

X 

4 

43 

O 

CP 

4) 

X 

4 

CM 

X 

4 

P- 

X 

4 

X 

X 

CP 

rH 

CM 

X 

4 

CP 

Ml 

fH 

p- 

X 

CP 

Ml 

fH 

<o 

4 

CO 

O' 

fH 

P^ 

4 

4 

fH 

4 

4 

X 

4 

43 

P- 

4 

CP 

O 

O 

X 

X 

CM 

P- 

rH 

X 

4 

CP 

X 

CM 

X 

4 

["> 

P- 

P- 

X 

X 

ir. 

4 

iC\ 

4 

X 

fyy 

4 

4 

X 

CO 

4 

X 

X 

P- 

X 

4 

tH 

X 

X 

X 

X 

CP 

X 

X 

CP 

X 

CM 

X 

f-l 

rH 

rH 

o 

O 

o 

o 

O 

o 

o 

O 

o 

O 

4 

4 

4 

4 

IT- 

4 

4 

X 

rH 

sO 

4 

4 

CM 

o 

43 

CP 

O' 

O 

X 

X 

X 

p'- 

p- 

P^ 

X 

P- 

P'- 

p- 

P- 

p- 

p* 

• 

p- 

• 

p- 

• 

• 

• 

4> 

• 

43 

• 

4) 

• 

4) 

• 

43 

• 

4) 

^o 

nO 

trs 

i£\ 

ITV 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

rH 

O 

c? 

o 

o 

CD 

c 

o 

o 

O 

rH 

fH 

fH 

tH 

rH 

fH 

t-H 

rH 

iH 

rH 

iH 

fH 

fH 

tH 

fH 

rH 

tH 

iH 

fH 

rH 

tH 

rH 

fH 

rH 

rH 

r-H 

rH 

tH 

rH 

o 

O 

o 

o 

o 

o 

o 

o 

o 

o 

o 

C' 

o 

o 

p 

o 

c- 

c 

c> 

Q 

O 

O 

O 

o 

O 

o 

o 

O 

O 

o 

o 

o 

o 

o 

o 

O 

o 

4' 

4^ 

4- 

4 

4 

4 

4 

4 

4 

1 

1 

1 

I 

1 

1 

1 

I 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

c 

o 

C- 

(D 

c. 

o 

o 

o 

o 

O 

o 

ID 

o 

o 

o 

o 

o 

Q 

(D 

o 

o 

o 

Q 

o 

o 

o 

o 

CD 

d 

CD 

Ci 

c 

o 

c. 

CD 

CD 

CD 

CD 

o 

o 

D 

o 

o 

o 

o 

o 

o 

X 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

O 

n 

o 

o 

c 

O 

X 

c 

o 

o 

o 

O 

-D 

O 

O 

o 

o 

o 

CD 

o 

o 

o 

o 

o 

o 

o 

o 

o 

CD 

o 

o 

o 

CD 

o 

o 

o 

CD 

O 

O 

Q 

o 

CD 

o 

o 

o 

o 

CD 

o 

o 

V 

CD 

CD 

o 

o 

O 

CD 

o 

o 

CD 

o 

CD 

o 

CD 

CD 

o 

o 

o 

<D 

o 

CD 

CD 

CD 

<D 

CD 

CD 

CD 

o 

V 

CD 

CD 

o 

o 

o 

o 

o 

CD 

o 

V 

v3 

CD 

CD 

o 

Q 

O 

o 

o 

o 

o 

o 

O 

o 

o 

o 

o 

o 

c 

o 

o 

o 

o 

O 

O 

O 

O 

O 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

V» 

O 

O 

o 

o 

o 

o 

O 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

O 

O 

O 

O 

O 

o 

o 

o 

o 

o 

o 

CD 

o 

o 

o 

o 

o 

o 

o 

CD 

o 

CD 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

Q 

o 

O 

O 

o 

CD 

o 

o 

'D 

o 

O 

o 

o 

p 

o 

Q 

o 

o 

O 

o 

c> 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

c 

o 

o 

o 

o 

o 

O 

o 

O 

O 

o 

o 

o 

o 

o 

o 

O 

O' 

X 

p'- 

X 

4 

<r 

X 

PJ 

tH 

o 

O 

o 

o 

o 

o 

CD 

o 

o 

o 

o 

o 

o 

CD 

CD 

CD 

o 

o 

Q 

CD 

O 

O 

o 

CD 

CP 

X 

p^ 

X 

4 

4 

U*! 

4 

4 

4 

4 

4 

4 

o 

o 

O 

o 

o 

o 

o 

CD 

o 

o 

o 

o 

o 

o 

o 

CP 

X 

P^ 

X 

4 

4 

X 

PJ 

rH 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

c 

O 

o 

o 

o 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

tH 

o 

X 

X 

4 

4 

X 

PM 

rH 

CP 

X 

p- 

X 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

66 


.OS'jCSOOO^'D-*!  4 .7044137a»0  +  00  1 .04399967b0-CA  4 . 1C3  18220  5D-0 1 


r 


•H 

*-4 

i-i 

•-I 

r4 

•-4 

f-4 

O 

o 

o 

O 

o 

o 

o 

o 

o 

• 

I 

1 

1 

1 

1 

1 

1 

1 

Q 

o 

e 

o 

o 

o 

o 

o 

o 

<M 

00 

<NJ 

lA 

lA 

Al 

GO 

lA 

•H 

m 

a:> 

CM 

CM 

A- 

o 

lA 

>o 

cx? 

o 

CM 

M- 

GO 

(X) 

CM 

>o 

f-4 

lA 

O' 

CO 

4> 

rH 

CD 

ro 

•H 

OO 

lA 

CO 

00 

O' 

CM 

lA 

GO 

o 

(A 

>o 

CO 

O 

O 

iH 

t-4 

CM 

CM 

CM 

CM 

o 

M 

o 

O 

O 

o 

O 

o 

r-l 

r-l 

1-4 

«H 

f-t 

<r 

>r 

■4' 

■J* 

M- 

M* 

M- 

43  A 

U>  A 

O 

<,» 

O 

O 

o 

o  o 

o  o 

o 

1 

1 

1 

1 

1  1 

1  1 

1 

O 

o 

o 

o 

o  o 

O  Q 

o 

CM 

a 

>r 

4> 

CO  A- 

<T»  03 

O' 

1—4 

M- 

M* 

^  A* 

U>  CO 

o 

CM 

<X- 

GO 

CM 

C5  CA 

M- 

o 

OO 

r-l 

CO  f- 

A-  ao 

M- 

fO 

A 

M- 

m 

00  ct 

A  ro 

•4* 

<r 

M* 

M- 

CA  M* 

no  CA 

OO 

M- 

tM 

sJ 

ro  fO 

<M  A- 

A- 

r— 1 

CO 

C^ 

o 

n 

43  >r 

CM  CA 

O 

o 

M- 

O 

CO 

W  CM 

CO  A 

cn 

r-4 

1 

A 

’  • 

• 

• 

• 

•  • 

•  • 

• 

o 

o 

o 

CM 

CM 

fH 

O'  U% 

r4  <M 

O' 

1 

O' 

i 

1 

1 

1 

1 

1  1 

1 

1 

o 

A 

Q 

o 

A 

O 

o 

CO 

o 

O' 

CM 

o 

A 

CM 

O 

o 

o 

O 

o 

O 

o  o 

C"  o 

o 

o 

f— 4 

o 

1 

f- 

O 

o 

o 

O 

o 

A 

• 

o 

o 

-4- 

4- 

-f 

4> 

4* 

4-  4^ 

+ 

O 

4- 

o 

o 

O 

O 

o 

o 

o  o 

O  O 

Q 

• 

« 

o 

CM 

CO 

•M 

r-4 

00  CO 

A  00 

A 

A 

rH 

o 

A 

A 

A 

>r  M- 

m  CM 

r-4 

o 

in 

m 

— 1 

O' 

r-  A 

ro  r-l 

H 

A 

o 

*D 

r- 

oc 

X> 

O'  r  > 

CM 

O 

o 

CM 

OC 

O  A- 

CO  O' 

CM 

M 

<. 

1  H 

o 

f-- 

^o 

>U 

O  A 

A 

A 

o 

O 

• 

M- 

M- 

M-  ^ 

sT  M* 

M" 

< 

CO 

4*  Z 

r4 

"7 

O 

O 

ri 

o  o 

fT  O 

CA 

1- 

x; 

4*  O 

A.  C- 

A- 

A. 

LU 

< 

^  —1 

• 

• 

• 

• 

•  • 

•  • 

• 

GO 

A 

>r 

M- 

M- 

M- 

'T  <r 

4-  M- 

4- 

A-  C/3 

CO 

U, 

UL 

O  Qu 

t>-l 

D 

o 

CT)  LU 

• 

UJ 

UJ 

IX; 

CA  Z 

A4 

1- 

I— 

•4 

f-H 

f-4 

►H 

-4 

r^  r-4 

rH  r-4 

rH 

< 

<I 

X 

C/3 

.1 

o 

O 

o 

o  o 

f  j  o 

O 

r: 

SC 

1— 

« 

1 

1 

1 

1  1 

1  1 

1 

*-» 

o 

O. 

Cl 

o 

o 

lH 

o  o 

o  o 

o 

1- 

A-  to 

LX 

> 

o 

r- 

O 

o  o 

rr  C3 

o 

t/O 

in 

(/3 

A- 

c»; 

Cj 

o 

o 

o  o 

o  o 

o 

IXi 

111 

LU  X 

to 

o 

o 

o 

o 

o  o 

o  o 

o 

O  -I 

o 

o 

o 

O  c 

w  Q 

o 

z 

«J 

■'5 

o 

c:> 

o 

o  o 

o  o 

CA 

< 

UJ  to 

■4 

OD 

r*- 

43 

M- 

CO  CM 

CM 

Z 

Z 

o  •— * 

Z 

O 

o 

O 

O 

o  o 

o 

O 

i—i 

OC 

f-H 

A 

A 

U3  A 

A  A 

A 

u. 

u. 

UJ  X 

u. 

o 

o 

o 

O 

O  CJ 

CJ  r> 

O 

>  CJ 

• 

• 

• 

• 

•  » 

•  • 

• 

UJ 

UJ 

z  *-• 

UJ 

A 

A 

A 

A 

A  A 

A  A 

A 

I 

X 

O  X 

X 

A- 

1— 

o 

►— 

67 


68 


PROGRAM  RUN  ENDS  NORMALLY 


Next  page  is  blank. 


69 


DISTRIBUTION  LIST 


Organization 

Commander 

Defense  Documentation  Center 
ATTN:  TCA 
Cameron  Station 
Alexandria,  VA  22314 

Commander 

US  Army  Materiel  Development 
and  Readiness  Command 
ATTN:  DRCCP 
5001  Eisenhower  Avenue 
Alexandria,  VA  22333 

Commander 

US  Army  Materiel  Development 
and  Readiness  Command 
ATTN:  DRCDE-F 
5001  Eisenhower  Avenue 
Alexandria,  VA  22333 

Commander 

US  Army  Materiel  Development 
and  Readiness  Command 
ATTN:  DRCRE-I 
5001  Eisenhower  Avenue 
Alexandria,  VA  22333 

Commander 

US  Army  Materiel  Development 
and  Readiness  Command 
ATTN:  DRCBSI-L 
5001  Eisenhower  Avenue 
Alexandria,  VA  22333 

Commander 

US  Army  Materiel  Development 
and  Readiness  Command 
ATTN:  DRCPA-S 
5001  Eisenhower  Avenue 
Alexandria,  VA  22333 


DISTRIBUTION  LIST 


No.  of 
Copies 

I 


1 


1 


1 


1 


1 
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Organization 

Commander 

US  Army  Materiel  Development 
and  Readiness  Command 
ATTN:  DRCQA 
5001  Eisenhower  Avenue 
Alexandria,  VA  22333 

Commander 

US  Army  Materiel  Development 
and  Readiness  Command 
ATTN:  DRCBSI-D 
5001  Eisenhower  Avenue 
Alexandria,  VA  22333 

Commander 

US  Army  Materiel  Development 
and  Readiness  Command 
ATTN:  DRCDE-R 
5001  Eisenhower  Avenue 
Alexandria,  VA  22333 

Commander 

US  Army  Materiel  Development 
and  Readiness  Command 
ATTN:  DRCDE-D 
5001  Eisenhower  Avenue 
Alexandria,  VA  22333 

Commander 

US  Army  Armament  Materiel  Readiness  Command 

ATTN:  DRSAR-SA 

Rock  Island,  IL  61299 

Commander 

Rock  Island  Arsenal 

ATTN:  Tech  Lib 

Rock  Island,  IL  61299 

Commander 

US  Army  Armament  Research  and 
Development  Command 
Dover,  NJ  07801 
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DISTRIBUTION  LIST 


Organization 

Commander 

Harry  Diamond  Laboratories 
ATTN:  DELHD-SAB 
2800  Powder  Mill  Road 
Adel  phi,  MD  20783 

Commander 

US  Army  Test  &  Evaluation  Command 

ATTN:  STEDP-MT-L 

Dugway  Proving  Ground,  MD  84022 

Commander 

US  Army  Aviation  Research 
and  Development  Command 
ATTN:  DRDAV-BC 
P.  0.  Box  209 
St.  Louis,  MO  63166 

Commander 

US  Army  Electronics  Research 
and  Development  Command 
ATTN:  DRDEL-SA 
Fort  Monmouth,  NJ  07703 

Director 

US  Army  TRADOC  Systems  Analysis  Activity 
ATTN:  ATAA-SL 
ATAA-T 

White  Sands  Missile  Range,  NM  88002 
Comnander 

US  Army  Missile  Command 

ATTN:  DRSMI-C 

Redstone  Arsenal ,  AL  35809 


Commander 

US  Army  Troop  Support  and  Aviation 
Materiel  Readiness  Command 
ATTN:  DRSTS-F 
4300  Goodfellow  Blvd 
St.  Louis,  MO  63120 
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DISTRIBUTION  LIST 


No.  of 

Copi es  Organization 

1  Commander 

US  Army  Tank-Automotive  Research 
and  Development  Command 
ATTN:  DRDTA-UL  (Tech  Lib) 

DRDTA-V 

Warren,  MI  48090 
1  Commander 

US  Army  Mobility  Equipment  Research 
and  Development  Coirmand 
ATTN:  DRDME-0 
Fort  Bel  voir,  VA  22060 

1  Commander 

US  Army  Natick  Research  and 
Development  Command 
ATTN:  DRDNA-0 
Natick,  MA  01760 

2  Chief 

Defense  Logistics  Studies 
Information  Exchange 
US  Army  Logistics  Management  Center 
ATTN:  DRXMC-D 
Fort  Lee,  VA  23801 

1  Commander 

US  Army  Concepts  Analysis  Agency 
8120  Woodmont  Avenue 
Bethesda,  MD  20014 

1  Commander 

US  Army  Logistics  Center 

ATTN:  ATCL-MER  (Mr.  H.  Burnette) 

Fort  Lee,  VA  23801 

Aberdeen  Proving  Ground 

Cdr,  USATECOM 

ATTN:  DRSTE,  DRSTE-CS-A 

Bldg  314 

Dir,  BRL,  ATTN:  STINFO  BR,  Bldg  305 
Director,  BRL,  Bldg  328 
Director,  HEL,  Bldg  520 
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